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Abstract” it is shown that the three types of diagrams discussed in Part | are members of a 


larger family, the defining feature of which is that they are homographically related, ice. trans 


formable into each other by central projection 


Phe rules of homographic transformation are outlined, and then used: (1) for the derivation 
of general rate-process constructions, (2) for the establishment of new types of diagram, the 
HO — air system serving as an example, (3) for the development of a new, entirely graphical 


construction procedure for determining the number of transfer units of exchange equipment 


\ new diagram is presented which is suitable for problems arising in drying practice, heating 


and ventilating and cooling-tower design An example of the latter is worked out 


Résumé — Loauteur a montré que les trois types de diagrammes ctudiés dans la part 1, font partic 
dune famille a caractéere homographique cest-a-dire transformable les uns dans les autres par 


une projection central 


Les révles de la transformation homographique sont esquissces et ensuite utilisces : (1) pour 


Veétablissement des constructions des processus gencraux de vitesse 


(2) pour Veétablissement de nouveaux types de diagrammes le svstéeme H,O-air servant 
d one 
(3) pour le développement un nouveau pr cede entiérement graphique de construction 


pour la determination du nombre des unites de trounsfert dun apparei’l dechange 


Un nouveau diagramme est presente qui est epproprie pour les probléemes se presentant dans 
la practique du séchage, du chauffage, et ce la ventilation et pour le calcul des tours de refroidis- 


sement. Un exemple de ce dernicr est étudic 


Zusammenfassung Die drei im Teil diskutierten Dingrammtypen sind Glieder einer grosseren 


Schar von Diagrammen. Das detinierende Kriterium ist die homographische Verwandtschaft, 


dh. die vegenseitige Transformierbarkeit durch Zentralprojektion, 


Die Regeln der nomographischen Transformation werden mitgeteilt und auf folgende Fille 


aunyvewandt 


1 Auf die Ableituny allvemeiner Konstruktionen fiir Ubertragungsvorgiinge. 

¥. Auf die Herstellung never Diagrammtypen, fiir die das System Wasser-Luft als Beispiel 
dient 

%. Auf die Entwicklung einer neuen, rein graphischen Konstruktion zur Bestimmung der 


Anzahl der Chertragungseinheiten cines Austauschers 


das besonders fiir 


Ferner wird ein neves Diagramm mitgeteilt Probleme der Trocknur 


Heizuny und Liftung und Kiihlturm-Berechnuny geeignet ist Fiir die letztere wird ein Beispiel 


ausyefiilrt 
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GENERALIZED CHARTS 
1.1 The Choice of Mass Unit 

We have seen, in Part I, [1] that graphical rate- 

process constructions can be carried out on the 

Hi — M chart just as well as on the A 

only their interpretations need modification, 

Now the H 


measuring quantity in terms of a unit (e.g. the 


m chart: 
M formulation is characterized by 


lb mole) w hich varies according to the « omposition 
of the mixture: the unit is 18 lb, for pure H,O, 


29 lb, for pure air, and has intermediate values 


for intermediate mixtures (see equation 6 of 
Part I) 


Examination of the arguments (not given here) 


Why should these units be chosen ? 


which lead to the construction of a chart on the 
‘ mole-of-mixture ~ basis, and to the formulae 
relating its co-ordinates to A and m. shows that 
at no point is any knowledge required of the 
actual relative masses of the molecules in question; 
even if matter were continuous (i.e. did not consist 
of molecules), the charts and constructions would 
remain valid, 

This recognition leads to an Important con- 
clusion: we may choose any size of mass unit 
whatever for each of the two components, Is 
there advantage in so doing’? Fig. 11, to be 
discussed later, shows a chart for H,O - air 
in which the H,O mass unit has been taken as 
1 lb,, and the air mass unit as 25 lb, (only the 
ratio of the units Is important, so one of them can 
conveniently be taken as unity). A glance at this 
new chart shows that the important region cover- 
ing gas-phase mixtures between 32°F and 212°F 


has been greatly enlarged; graphical construc- 
tions for these mixtures are thus easier to make. 
So the freedom to choose the mass-unit ratio 
makes it possible to * blow-up particular 
regions, while still retaining the whole range of 
mixture compositions (pure air to pure H,O) 


on the chart. 


Mollier 


“ Molecular weight” of water on the 
(1 — x) chart 


It is interesting to note that Mouiier has 
implicitly made extreme use of this freedom in 
in constructing the 7] —v chart: he has made the 


“molecular weight” of water, tend to 
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infinity. kor in this case, equation (5) of Part I, 
for example, takes the form (with ay 4 in the 
place of 18, and 1 in the place of 29); , 


M m (1 m) (1) 
Comparison of (1) above with (5) of Part I shows 
that M is the same asx. So the J chart 
is really an H-—M chart, but with Mee = 1, 


= Very large. (The appearance of On 
the left-hand side of (1) is insignificant, for we are 
always free to introduce a constant multiplier 
when choosing a seale to which to plot the 
diagram). This choice of mass unit happens sO 
to cramp the region of moderate H,O contents 
that the part of the diagram dealing with large 


HO contents has in practice to be discarded. 


Further generalization 


Although the free choice of mass unit already 
provides possibilities of chart design which will 
not be quickly exhausted, an even wider range is 
possible. These are revealed by considering the 
charts from the point of view of projective 
geometry. We first state the features of this 


branch of mathematics which will be used below. 


1.2. Some Facts about Projective Geometry 

(i) Projective geometry concerns the relations 
between a diagram drawn on one plane and its 
shadow on a second plane when illuminated by a 
point source of light. The process of forming the 
second diagram from the first is known as central 
projection. 

(ii) If the co-ordinates of a point on the first 
plane are (/, m) in a rectangular system, and the 
co-ordinates of the corresponding point on the 
second plane are (», €) in a second rectangular 


system, the co-ordinates are related by 

a,m m bah (2) 
pm gh +r pm gh +r 


the relative positions of the two planes and the 
light source (vertex of projection), and on the 
chosen co-ordinate axes, and so are valid for all 
parts of corresponding points on the two planes. 
Note that the denominators of the two expressions 
are identical. 
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If (2) is treated as representing an analytical 
operation, it is known as homographic transforma- 
tion, 
(iil) 
in the h 


in the »-€ plane, and vice versa. 


In a homographic transformation, a point 

m plane corresponds to just one point 
Similarly, 
there is a one-to-one correspondence between 
straight lines on the two planes. 

(iv) If four points A, B, C, D lie on a straight 
line (not necessarily in the order stated), the 
product of the length ratios 


AB. CD 
AD. CB 
is known as the cross-ratio of the range of points 
A, B, C, D. 
(v) If the straight line A, B, C, D is trans- 


B’, C’, D’ homo- 
by central projection, the cross- 


formed into the line straight A’, 
graphically, i.e. 
ratios of the two corresponding ranges are equal. 

(vi) Ifa range of points A, B.... P,Q.... 
on a straight line / is in homographic relation to 
second straight line /', then the intersection of the 
lines PQ’ and P’ Q lies on a fixed straight line 
which passes through the point of intersection, 
A, of l and l' (Fig. 1). 

(vil) 


homographic relation, in 


a range of points A, B’,.... on a 


If two planes h-m and »-€ are in 
general one of the 
straight lines on the —m plane does not appear 
‘vanishes or 


on the »- € plane: it 


This is the 


goes to 


infinity.” line pm +gh +r = 0. 


re) 

Pr 

B P g 


Fic. 1. Ulustrating that the intersection of PQ’ and 


PQ lies on a fixed straight line. 


Conserved-property diagrams for rate-process calculations 


Similarly one of the ines on the »~— & plane 
corresponds to the line at infinity on the h—m 
plane. 
(viii) 
on the h 
€—» plane and vice versa, 
(1x) 


Conics (ellipses, parabolae, hyperbolae) 


m plane project into conics on the 


Any conic can be projected into a circle 
having the projection of any given point (not on 
the conic) as its centre. 

Of these statements, the most important for 
the present purposes are (ii), (v) and (vii). All the 
statements are used explicitly or implicitly, in 
the following pages, 


1.3 Application of Projective Geometry to the 
Conserved-Property Charts 


Relevance to the h-m, H- M, and 1- charts 
for H,O — air 
Comparison of equations (4), (5), (10) and 


(11) of Part I, 
homographic transformation (2) above, shows that 
the h—-m, and 
graphic relation to each other, with M and w& as 
particular forms of € and H and TJ as particular 
forms of ». 


for example, with the general 


charts are in homo- 


The transformations relating the three charts 
are not of the most general form since only the 
horizontal co-ordinate appears in the denominator 
of the equations. It is for this reason that lines 
of constant composition remain vertical in each 
of the charts. The transformations are further 
restricted in that M 0 when m 0. 

In the transformation between the h 
H—-M only the co-ordinate 
appears in the numerator of the transforming 
relation (10) of Part I; so the line corresponding 
to h both charts. In the 
transformation from A to J on the other hand. 


0 and a 
m and 


planes, vertical 


0 is horizontal in 


the horizontal co-ordinate m does appear in the 
numerator (equation 11 of Part I); so the line 
h 0 (or 7 0) becomes a sloping line on the 
I—x chart. This has already been commented 
on above. 
The line m = 1 (pure H,O) on the h—m chart 
is the line at infinity on the J— 2 chart. In the 
from the h—m to the H—~M 
plane, it is the line (18 29) + [1 — (18/29)] 


m 0 does 


transformation 


which vanishes: since this not 
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correspond to a realistic composition (for m cannot 
in practice be negative), no physically interesting 
mixtures are lost. 

In general, the line which goes to infinity in a 
transformation to the H—M that 
* molecular 


plane Is 
corresponding to mixtures of infinite 


weight” p. 


The possibility of other enthalpy-com pos fion charts 


The chief 
enthalpy-composition charts can be generalized, 


wav in which the construction of 


in addition to the free choice of already 


mentioned. lies in introducing the vertical co 


ordinate h into the denominator of the transform- 
This 


composition ceasing tobe vertical, or even parallel ; 


ing relations. results in lines of constant 


the region of realistic mixtures maps, in gr neral, 
on to a triangular area. Interpretation in terms 
of the choice of mass unit becomes more diflicult : 
enthalpy is now being treated in some rv spects as 
equivalent to mass. 

There are nine constants in the general trans- 
formation (2) [corresponding to the six degrees 


of freedom in choosing the position ot the second 


plane relative to the first plus the three degrees of 


freedom in choosing the position of the vertex 
of projection (light-source)]. However, only four 
of these are of great interest ; the other five merely 
govern the horizontal and vertical positions, the 
horizontal and vertical scales and the orientation 
of the axes of the graph paper on which the 
transformed chart is being drawn. Of the four 
interesting constants, the two in the denominator 
of the expressions for » and € determin which 
the h—m chart 


infinity ; the two remaining ones determine the 


line on will bn projected to 


slopes of the projections on the plane of the 


m 0 and h 0 lines. 


Summarizing, the enthalpy-composition chart 
for mixtures of a given pair of substances at a 


given pressure can be drawn in a four-fold 


infinity of ways, mere scale changes being left 
out of account. This allowance should be suflicient 


for most purposes 


present pair of papers the 
nized the 


*Since writing the 
This work shows clearly that BusEMANN rece 
of diagram, although he 


heation appears to have been uncese rvedly me u lec ted 


author has obtained the 
projective 


did not exploit his recognition in the 
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Diagrams for ternary miatures 


Diagrams showing the phase boundaries and 
mixed-phase isotherms of ternary mixtures are 
similar in many respects to enthalpy composition 
diagrams for binary mixtures: the mass or mole 
fraction of the third component takes the place 
of enthalpy. The chemical engineering literature 
the 


three compone nt 


triangular diagrams (in which 
that of the 


solvent-free basis *’) 


abounds in 
quantity unit ts 
mixture), and rectangular ( 
diagrams (in which the quantity unit ts that of 
only two of three components). 

It should be 
that the triangular diagrams and the solvent-free 
that 


exists in the mass units used, 


clear. without detailed discussion, 


ones are homographically related, vreat 
freedom of choice 
that the triangular diagrams do not have to be 
equilateral triangles, and so on. Ternary-mixture 
exhibit special features 


diagrams do not any 


therefore. 


The Lever Ruk 

The 
property charts is the Lever Rule, 
of Part 1. 


expressed as a cross-ratio ; 


in terms of cross-ratios 


foundation of the usefulness of conserved- 
equation (1) 
It will now be shown how this can be 
the result will then be 
used to provide simple proofs of some of the rate- 
process formulae. 
Consider the h 
X and Z represent systems which are mixed to 
form a new system Y; F 


» Z with the 


m chart of Fig. 2(a). The points 


is the intersection of 


line m 0. From equation (12) 


of Part 1, and consideration of Fig. 2(a), we can 
write : 
Wz XY.ZF omy 
wx YZ XF.ZY mz 
W7 XY . ZF (3) 
Wx MX ZY 


Now the expression on the right-hand side of 
the 


cross-ratios 


this equation is the cross-ratio of range of 


points X, Y, Z, F. So, 


under 


since are 


transformation, 


invariant homographic 


book by BusemMann referred to in Part 1. 


of the connexion between the various types 


nature 


manner of the present work BuseMANN's pub- 


: 

l] 

‘ 

950 

ake 

= 
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m=O 
h 
Fe. x 
~~. 
— 
| Z 
en 
| 
(a) (b) 
hic. 2. The Lever Rule on the Y chart 
X¥ X\¥ wz om (éz ) 
VOL, (a) (b) : 
Wy Z Wx my Z 
os f) 
the expression may be evaluated on the n-€& validity of the procedure described under Section 
chart (Fig. 2b) formed from the h— m chart by 2.4 (b) (iii) of Part I, 

. central projection. Now on the » - € chart, the In calculating the change of gas state resulting 
ratio ZF/XF is equal to (fy - fp) (Ex — &p). from a contact between gas and liquid, it was 
Equation (3) therefore can be written: shown above that the increment in the number of 

—_= gas-side mass transfer units was equal to the 
W7 (Ex Cp) = = . 
| = ay (4) quantity G,G,.TJ/G,J.TG,, (equation 21, 
YZ wx Mx lee vid Fig. 5 of Part I). This quantity may now be 

4 in which the suflix »— & signifies evaluation on recognized as a cross-ratio: it therefore does not 

the » — € chart. matter on which diagram the lengths are 

If now we restrict attention to transformations evaluated. This was why the same procedure 
of the HM variety (which includes 7—2) in could be used on the H—M chart as on the 
which the horizontal co-ordinate (i.e. M) is chart. 


equal to zero whenever m 0. equation (4) 


reduces to equation (13) of Part I, 


Rate-process constructions in terms of cross-ratios 


The constructions involved in heat and mass 
transfer rate calculations require the evaluation 
of length ratios for which one point, for example 
H in the ratio AG ‘HA, lies on the line m s 
We how AG HA on an H—~M 


chart is related to that for the corresponding 


will now see 
h—m chart. 
Appropriate substitutions — in (3) 


(M for & H for », M, 


and (4), 
0, G for X, H for Y. 


A for Z, my My 1), lead to 
AG AG M 
(5) 
HA/ Me 


which, taking account of (6) of Part I, proves the 


Part 


Miztures 


-phase 


14 Some 


with 


Charts for H,¢ )— Air 
parallel 


Diagram horizontal 


isotherms 

The specific heats at constant pressure of steam 
and air at moderate temperatures are in the 
0-445 : O24. As a of this, 


and of the choice of enthalpy base for the h — m 


ratio consequence 


chart, all the gas-phase isotherms converge on 


the point (m*, h*) where m* 1:17 and 
h* 1260 B.t.u ‘Ib,,. 

In order to obtain a diagram with parallel 
gas-phase isotherms, we project the line m m* 


to infinity. This is effected by the transformation : 


*The symbols 7 and M are used instead of » and & 
wherever the constant-composition lines are vertical and 
the whole range of compositions appear in a finite hori 


zontal range 
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2-17 m 
Mt 6) 
1-17 


In order to correspond with (6), and in addition 


to cause the gas-phas« isotherms to be horizontal, 


the following transformation of the vertical co- 


ordinate can be used 


2-17 (h 


mh, 
a 


m 1-17 


H (7) 

The corresponding diagram is shown i Fig. 3. 
It has the advantage that the gas-phase part of 
the diagram is particularly easy to draw. 
Physically the 
thought of as employing mass units of 1-85 Ib. 


of air and 1 Ib,, of H,O. 


interpreted, diagram may be 


VW chart for HO irat ate 


parable was hase teothertius 


pre sslire 


Diagram using mass of H.O as basis 


The inability of the I ~x diagram to represent 


verv moist mixtures has been noted. If condensa- 


tion of steam containing small quantities of air 
is to be studied, it may be prefe rable to project 
the m Th 
transformation of the horizontal co-ordinate ts: 


0 line to imfinity. corresponding 


(1 m)m (3) 


A suitable transformation for the vertical co- 
ordinate 1s 
(9) 


200°F / 


Woter line 


50 


E.nthalpy-composition chart for H,O - air at 


latm pressure basis. 


This makes the intersections of the isotherms with 
the 


line. 


0 line identical with those on the m l 


The resulting diagram is shown in Fig. 4. Of 
course only a limited range of € can be used ; 
the range 0 to 5-0 has been chosen in the example. 
The diagram may be interpreted as employing a 
mass unit for air which is very much larger than 
that for H,O. 

N.B 
at (— 2171-17 
drawing the diagram) 


The focus of the gas phase isotherms (£*, »*) is now 


1260 1-17); plotting this point assists in 


A triangular diagram 


To underline the fact that enthalpy-composition 
and ternary-mixture diagrams are all members of 
the same family, a transformation will be intro- 
duced which maps the h—m plane for H,O-air 
on a triangular area, which may be equilateral. 
A transformation giving an isosceles triangle 
is: 

bh h/hy, (10) 
1 + hyhy, h 


hy, can be taken as 1076 B.t.u/Ib,, H,O. 
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\ triangular enthalpy-composition diagram for H,O — air at 


Fig. 5 shows the resulting diagram ; the vertical 
and horizontal scales have been chosen to make 
the vertex angle of the triangle 60°. Since points 
of infinite enthalpy appear on the diagram (at the 
vertex), such charts may be useful when very 
dealt 


with, for example those obtaining when a body 


large temperature differences have to be 


travelling at high Mach Number is cooled by 
Note that the focus 


of the gas-phase isotherms is at £* ‘BO, 


water flowing on its surface. 


” 6-86. This is far off the diagram; so 


the gas-phase isotherms are nearly parallel. 


Other diagrams 


As will by now be apparent, the number of 


conserved property diagrams which can be drawn 
is inexhaustible: one must choose the diagram 
which has the most convenient shape for the 
particular computational problem = which has 
arisen. The following suggestions may be worth 


making here however. 


diagrams for rate Process « theulations 


Part 


atm pressure 


Mixtures with large moisture contents do not 
have to be represented on diagrams such as that 
of Fig. 4; this diagram, like the J —a 
the disadvantage of lacking one 
Instead it may be preferable to use an H—M 
transformation employing, in effect, a very small 


one, has 


boundary 


mass unit for air compared with that for H,O. 
This has the effect of stretching the part of the 
h—m diagram near m 1 and contracting the 
part near m = 0. 

If a suitably small mass-unit ratio is chosen, 
it is possible to exhibit, as an area of finite size, 
the region of air-in-water solutions. The diagram 
can then be used for solving problems connected 
with de-aerating plant, and the like. 

The plotting of the saturated-gas line is the most 
tedious part of the construction of an enthalpy- 
composition chart. Yet if the interesting part 
of this line can be approximated on the h-m 


plane by a part of a conic, projective geometry 


shows that it is possible to find a homographic 
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transformation which projects this come into a 
circle with any desired point as centre, so pe rmitt- 
ing the saturation line to be drawn with compasses, 
the 
this should permit the devising of 
the 


Further possibilities lie in Principle of 
Duality [2] 
diagrams embodying thermodynamic pro 
perties of a mixture system, amenable to con 
structions such as those of the present paper, 
in which the isotherms are repres nted by points 
Whether the resulting simph- 
the 


of such a representation can only be decided by 


rather than lines. 


fication can counter-balance unfamiliarity 


further study. 
It is troublesome to have to plot a new diagram 
therefore 


for each pair of substances. It seems 


desirable to study whether, by suitable choice of 


reduced co-ordinates, the equilibrium diagrams 
for various pairs of substances can be made to 


The 


of co-ordinate presented by projective geometry 


correspond closely. new freedom ol choice 


may make this possible. 


RatrE-Process CONSTRUCTIONS ON THI 


GENERAL CoNSERVED-Property CHAR’ 


The constructions will now be given which 


caleu- 


pe rmit heat and mass transfer rates to be 


lated with the aid of » charts. These may be 


Fig. 6 Rate 


revarded as generalizations of those in Section 
2.4 of Part I. The H,O air system will continue 
to be used as an example, but no restriction of 
generality is implied thereby except, where 
straightness of gas-phase isotherms is assumed, 
to ideal mixtures. 

Thereafter, new constructions will be presented 
which permit graphical determination of the size 
of equipment needed to perform a given transfer 


operation. 


21 Determination of the Interface States Land S$ 
(Fig. 6) 
The problem and notation are as for Section 2.4 
above. The procedure is : 
(i) Locate O on the pure HO line corre spond 
ing to water at f,. 
(ii) Locate H on the pure-HH,C ) line correspond- 


ing to water at f,,, which is given by 
(11) 


(iii) Draw the line HG and tind A on it such 


that 


process and gas-state change construction. 
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where F lies on HG and on the pure-air (m 0) 
line, 

(iv) Locate L and S which lie at the extre- 
mities of the mixed-phase isotherm through A. 
L, lies on the pure-H,O line ; 8 lies on the satura- 


tion line. 


Transfer Rates 


( onstruction 


2.2 Determination of (Fig. 6) 


(i) Loeate D on the isotherm 


through G such that may Me. 


Cus phase 


N.B. Lines of constant m are conveniently drawn on the 
chart Of course, when &€ depends only on m these lines 


ure vertical in triangular diagrams like Fig. 5, they all 


pass through the vertex of the triangle 


(ii) Produce LS to cut DG at C. Plot K at the 


intersection of DG and the gas-phase isotherm 
tx. where 


(13) 


Transfer rates 


The mass transfer rate is given by 


( s ye ao ( &¢ fy 
(14) 
SI, mM (oy fy ) 
where F’ lies on CL and on the pure air line 
The heat transfer rate to the reservoir is given 
by 
a 
hy h, (15) 


The heat transfer rate from the gas to the 


liquid surface is given by 


h hy (16) 


N.B. For the last two evaluations, it is convenient to 


have lines of constant A drawn on the chart 


2.3 Change of Gas Condition (Fig. 6) 
Construction 
(1) Produce KS to cut the pure-H,0 line at J. 
(1) Locate T where G,J cuts the constant- 
composition line through S. 


N.B. G is now called 
the contact surface. 


it represents gas approaching 
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(ii) Find G, on G,J such that 
GG, . TJ 
GJ. TG, 

Comments 
Sometimes it is preferable to fix G,, arbitrarily, 
at, say, some fraction of the distance from 
G, to T. 


the left-hand side of (17) gives the necessary 


Then evaluation of the expression on 


contact area in terms of number of gas-side mass 


transfer units. BoSnsgakovic [3], for example, 
in his constructions carried out on the J —a2 
chart, suggests making G,G, G,T 1 2 for rough 
calculations. 

Kquation (17) is, in essence, a finite-difference 
approximation to a differential equation. It is 


when is 


therefore only strictly valid 


infinitesimal. Except in very rare circumstances 
however the finite-difference form has to be used 
in order to permit integration, i.e. determination 
of the number of transfer units of the equipment. 
Then the usual compromise between truncation 
and rounding-off errors has to be made in choosing 


the size of step. 


2.4 Determining Corresponding Gas and Liquid 
States in Transfer Equipme nt 

So far we have considered heat and mass 
transfer at an infinitesimal element of the contact 
surface. In plant such as cooling-towers or packed 
distillation-columns, both gas and liquid states 
vary with position; the corresponding interface 
states vary also. We now turn to the graphical 
procedures for determining these sets of states, 

Part of the technique for doing so is well 
known: it involves the use of the so-called polar 
diagram [4, 5, 6]. For adiabatic parallel-flow or 
counter-flow plant, this diagram consists of a 
pencil of lines radiating from a single point on a 
conserved-property diagram, the position of 
which is fixed by the states of the ingoing fluid 
streams and the ratio of their mass flow rates. 
The importance of the polar diagram lies in the 
fact that the state-points for the liquid and gas 
streams, at any level of the equipment, lie on 
one of these lines. The validity of the construc- 
tion derives from the laws of conservation of mass 
and energy. 


Satn line 
y 


(a) for adiabatic equipment 


Fig. 7(a) shows such a polar diagram for the 
adiabatic counter-flow plant shown in Fig. 7(b). 
G represents the gas-stream and W the liquid. 
G,W,, G,W, ete. all pass through the ~ 
P., which represents the state of a 


counter- 
flow pole,” 
single stream from which the pairs of liquid and 
gas streams could be derived. 


The position of P is such that 
me, 


PG, 


Moy (fw, n 


(18) 


(SGn 


where the subseript n denotes the tower section 
in question, mi, and ry denote respectively the 
mass flow rates of the gas and liquid streams, 
and F,, is the intersection of the line P_W,G, 
with the pure-air line. 

Normally the pole P_ is determined from (18) 
through specification of the liquid and gas states 
at one end of the plant, say, W, and G,, and of 
their mass flow rate ratio there. Then the rays 
P.W,, P.W,, ete. are drawn. 

In order to locate the points G,, G,, etc., it is 
common to assume that they lie on the saturation 
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line ; they are then completely fixed by the rays 
through P.. In a rigorous calculation, however, 
it is necessary to apply the construction for the 
direction of change of the gas state several times 
in succession, starting at G, and finishing (in 
the example shown) at G,. The 
for the gas in the equipment, 


result is an 
‘operating line 
i.e. a curve which gives the gas-state for any 
water state. 

The analysis of the equipment thus proceeds 
in three stages : 

(i) P. is fixed. Then conservation principles 
ensure that the G-states lie on the appropriate 
PW lines. 

(ii) Rate-process considerations, giving the 
direction of movement of G, combine with (i) to 
allow the position of the G’s to be fixed. 

(iii) 


to determination of the number of transfer units 


Further rate-process considerations lead 


of the equipment. 
Suflicient attention has already been given to 


stages (i) and (ii). We now turn to stage (iii). 


= 
| 
. 
F —_ | 
Ad ~*~ G : 
Ww 
3 
3 mi | V 
| 
’ 
hie. 7. 
7. 
a 
5 
|| 
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(N.T.U.) 


2.5 Determining the Equipment Siz 
The problem 

The column height of the transfer process shown 
in Fig. 7(a) can be obtained simply by evaluating 
the expression (17) for each of the stages shown. 
This procedure is straightforward but involves 
the 
ratios. 


measuring of lengths and 
It would be more convenient if a purely 
graphical integration procedure could be devised, 
so that the only numerical step occurs at the very 


end. 


Graphical integration procedure 
The 


described. Then its validity will be demonstrated, 


recommended procedure will first be 


and other possibilities discussed. 


Fig. 8 shows a part of the construction for five 


gas states. G,, G,...... G 


. represent succes- 
sive gas states, and the corresponding points T,, 
those featuring in the rate-process construction 


are 


and in equation (17). The steps are as follows: 
(i) Draw any line cutting the pure-H,O line 


at j. 


(ii) Draw parallels to the pure-H,O line 


through the G’s and T’s, to intersect the line just 
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Part Il 


m,2° 


uit gst)] 


drawn in G,’, G,’...G,’.... and T’,, T,’..... 

(iii) Plot T,,’ at the mid-point of T,’T,’; 
plot T,,' at the mid-point of T,’ T, ; ete. 


(iv) 
Choose any point, g,, on it. 


| 


Draw a second line through the point j. 
Plot t at the mid- 
point of g, j. 

(v) Project the range of points G’,, G’,, T's 3 
into the range of points g,, g», t, j. This fixes the 
point g 


(N.B. 
simpler is shown in Fig. 9 (a): 


The 
and are 


Two methods of projection can be used. 
g, G, 
produced to cut at R ; RG,’ is produced to cut g, j at 
This construction is inconvenient when g,G,’ and 
t Tye are nearly parallel ; for then R may lie out of reach. 
An alternative construction, which does not suffer from 
this defect. is shown in Fig. 9 (b): 
drawn N; g, Gg 

jN at Q; G,’ Q is drawn to cut g, j at go) 


and g Ty are 


| 


intersecting at is drawn, intersecting 


(vi) Project the range of points G,’, G,’, T,,’, j 
into the range gp, gs, t, j. This fixes the point gy. 
(vii) 


tions of the same procedure. 


Find any other g,, by successive applica- 


(viii) Measure the ratio of lengths g, t/g, t. 
(ix) Evaluate the number of gas-side mass 


. J 
rae, 
* 
G, | 
G 
VOL, 


D. B 


t 


Fig 9. Two 


transfer units necessary between level 1 and level 


n from 


N Inj} (1 g,t (19) 


(NB 


have been drawn on a suitably 


The last two steps are avoided if the points 


marked scule vided on 


the diagram. Only one such scale is needed 


otherwise free to choose the point j anywhere on the 


pure — H,O line) 


Validity of the construction 


The range of points G,’, G,’, T,’, } is obtained 
from the range G,, G,, T,, J, by central projection 
(the vertex of projection happens to be at infinity 
in the recommended procedure, but it need not 
be); therefore the cross-ratios of these ranges are 
The 

From 


is true of the other 
the 


is given by 


same pairs ot 


(17) 


gas-side mass transfer units N 


equal, 


ranges. therefore, number of 


m.2 


G 


n 


(V0) 
G 
The sutlix for T in this equation could be either 


n or n 1; for the equation is only strictly 
valid when the difference between the n th and 
(n + 1) th With ste ps of 


finite size, it seems wisest to use the corre spond. 


states is very small. 


ing to the mid-point of the interval, namely 
The 
from the range G,, G,’, T,,’, } by central projec- 
Their 
The same is true of the other pairs of ranges. 
From (20), therefore, 


range of points g,, t, is derived 


tion. cross-ratios are therefor equal, 


(21) 
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ethods of determining the position of g 


Now for a very large number of intervals, (21) 


can be written 


This 
since & and &, 
has been specially designed to tix these points). 


We obtain 


integral can be evaluated analytically, 


are constants (the construction 


(23) 


For the particular case in which t is midway 


between g above, this reduces to 


| 
equation (19), 
It is that the 


procedure gives a finite-difference approximation 


and j, as 


therefore clear recommended 


to the solution of a differential equation; it gives 
an exact solution, if drawing errors are neglected, 


when the size of step is small. 


Comments on the procedure 


Many 


be devised, 


the construction can be 


For example, t need not be at the 


variants of 


mid-point of g,). though this arrangement gives 


the greatest accuracy if g, } 1s finite. N can be 


read from the scale with still higher accuracy if 
} is projected to infinity, though this slightly 
complicates the construction. The first projection 
(G,, to G,’) need not be a parallel one ; any point 
on the pure-H,O line can be taken as vertex. 
If all the T’s lie approximately on a straight line, 


there is advantage in making the vertex just 


mentioned lie on this line also. Many other 
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: 


possibilities will occur to anyone performing such 
constructions. 

The time taken to carry out an integration, 
once the G's and T's have been determined. 
does not exceed more than a few minutes in most 
practical problems (up to, say, ten steps). The 
mental effort, after some practice, is also small. 

The accuracy of the result falls off as N , 
increases and will often be unace« ptably low for 


greater than 3. nt involving large 
numbers of transfer units is probably best cal- 
culated in two or mor parts, the total size 
being caleulated by adding the N 


An example of the use of the recommended 


m 5 of the parts. 
procedure is given below in connexion with a 


cooling-tower. 


3. A Useru, Cuarr ror H,O- Air Mixtures 
3.1 Th scription of the Chart 


For many of the problems arising in heating 
and ventilating, industrial drving, and cooling 
tower design, the ** mass-of-mixture ” basis \ ields 
enthalpy-composition charts on which the state- 
points of interest are restricted to an inconven- 
iently small area. The chart shown in Fig. 10 is 


designed to avoid this difliculty,. 


The transformation 

The chart is of the H — M variety, i.e. lines of 
constant composition remain vertical, and the 
horizontal range is from 0 to 1. The co-ordinate 


transformations are 


M - - (24) 
(1/25) (24 25) n 
h 
H - (25) 
(1 25) (24 25) m 


The numbers appearing in the denominator are 
chosen so as to make the gas-phase isotherm 
approximately horizontal at the atmospheric 
boiling-point (212°F). It will be noted that, 


close to m l 


M~m; li~h (26) 


close to m 0 


M~25m; H~25h (27) 
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So a 25-fold enlargement of the air side of the 
diagram has been effected. when compared with 
the h — m chart, while the H,O side of the diagram 
is unaffected. 


The MASS has is 


The transformations (24) and (25) permit 
interpretation as involving the use of 1 lb, as the 
mass unit for H,O but 25 lb, as the mass unit for 
air. The effective ** molecular weight ” is thus 

I 24 
p=1/(— m (28) 
J 


23 


For the sake of simplicity of speech, , Ib, of 
mixture is called 1 **lb nol”’ on the ordinate scale. 
Seales of «, and of m, are provided at the base 
of the diagram. .M (and m) happen to increase 
from right to left on Fig. 10, to accord with the 
practice of [7] and [8]; this fact is without other 


significance. 


Saturation lines 


Provided that air and steam obey the Gibbs- 
Dalton laws, and that, like water. their enthalpies 
depend on temperature alone, which may be 
assumed for practical purposes, the positions of 
the gas-phase isotherms are independent of 
pressure, It is convenient therefore to make a 
single chart valid for several pressures by drawing 
several saturation lines. Those for 0-05, 0-1, 0-2. 
0-4, 0-6, 1:0, 2-0, 5-0 and 10-0 atmospheres are 
drawn on Fig. 10. 

To use the chart for mixtures at a particular 
pressure, the appropriate saturation line is 
chosen. The branches of the gas-phase isotherms 
to the left of their intersection with the saturation 
line are now ignored. In their place are drawn 
mixed-phase isotherms ; these are straight lines, 
terminating at one end on the temperature scale 
marked on the water line (M 1), and at the 
other end at the appropriate intersection of the 
gas-phase isotherm and the saturation line. 

The saturation line for 1 atm (14-696 Ib, in?) 
total pressure is distinguished by hatching. If its 
shape is compared with that on an h—m chart 
(Fig. 2a of Part I), it will be seen to be more 
curved; a point of any given temperature is 
displaced towards the line m 1. 
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Enathalpy-composition 
1 nol air 


Fic. 10 
nol H,O Ib, H,O 
= enthalpy nol mixture m 


(1 25) 


Numbers underlined at lower end of saturation lines represent 


Number-of-transfer-units scale 

On the base of the diagram is a scale for use in 
calculations, plotted in accordance with 
equation (2: 1, g, on the 
line M = 0, and t on the line 2. The 
numbers give the value of N,,. when the final 
gas-phase point g, lands on the corresponding 


Nas 
8); j is on the line M 


graduation mark. 


25 Ib, air 


24 25 


M) chart for H,O-air at various pressures. 


M nol H,O nol mixture 


il 
total pressures in atm. 


(24 25) 25) 24,25) 


Applications of the chart 
The extension of the low-m part of the diagram 


is valuable whenever the water content of air at 
in question, as in 
For example the 
the 


moderate temperatures is 
psychrometry and in drying. 
composition of the air entering 
tower in the next example is easily fixed, knowing 


the wet-bulb and dry-bulb temperatures. 


cooling- 
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3.2. A Cooling-tower Example Data 
To illustrate the use of the chart, and of the eal- Water inlet temperature: 110°F 


culation method of section 2.5 of this paper, we pre- Water outlet temperature: 85°F 


sent an example of the determination of the number Air inlet condition ; dry-bulb temperature : 8: 
of transfer units (N.T.U.) of a cooling-tower. wet-bulb temperature : 
The complete construction is shown in Fig. 11. Total pressure: 1 atm 
Air/water ratio: 
1-5 times the minimum. 


Operating ine 


B.tu. nol 


4, 


G2 


Fic. 11. Determination of number of gas-side mass transfer units for cooling-tower example. 
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B. SPALDING 


and Las 


Determination of corresponding water 
stream states 
the 


are 


The state-point of the water at inlet W,, 
water at outlet W, and the air at inlet, G,, 
marked on the chart (Fig. 11) directly from the 


data. 


N.B. The suffix 6 to W at inlet is a consequence of 


deciding, in addition, to consider conditions at six levels 


in the tower. characterized by 5 F differences of water 


temperature 


The minimum permissible ratio of air rate to 
water rate is that causing the air to leave the tower 


in equilibrium with the incoming water: so 
110 F. 


G is on the l atm saturation line at 


6 min 
The 
is then 


and W,G,. 


corre sponding 
drawn as the imtersection of WG 


We have 


6 thn 


mn 


Mey 


Me, 


(20) 
p min 


Since the air rate is to be 1-5 times the minimum, 


we find P,. such that 


W, P, 
P, G, 


l 5 6 


(30) 
P, G, 
This is done by means of dividers. P,, is now the 
parallel-tlow pole of the tower. 

W, 


intersection is marked G, (for we know that the 


is produced to the saturation line; its 


air will be practically saturated when it leaves 
the tower). W, G, and W, G, are now drawn to 
intersect at P.. 
of the tower; it has the property that the state- 


This is the counter-tlow pole 


points for the water and air streams at any level 
in the tower lie on lines passing through P.. 
PW, Gy, 
PW, G, P, W, G, and P, W, The 
points W,, W,, W, and W, are on the liquid line 
at 90°F, 95 F, 100°F 05°F. The 
G,, Gs, Gy, G, cannot yet be finally located, but 


Four more such lines are drawn : 


and points 


we know that they lie on these lines. 


(N.B 


actually less accurate when carried out on Fig 


The part of the construction just des ribed is 
11 than on 
diagram, or an // —- M diagram 


1 for the 


a straightforward 


with a smaller mass-unit ratio than 25 points 
r, and i. on Fig. 11 turn out to be rather near to the 
W's. Indeed there 


for the whole of the N.T.1 


is nothing wrong with an  m chart 


construction, proy ided that the 


vertical scale is suitable Hlowever it is easy to derive 


equivalent: constructions which retain high accuracy on 


Fig. 11 


us well) 


Las ope rating line 


Determination of the 


It will be 
tower design, that (1) o i. 


assumed, as is common In cooling- 
(ii) the liquid-side 
The 
consequences are that the gas-phase surface states 


5S, 


6 


heat-transfer resistance is very small. 
can immediately be plotted at the 
other extremities of the mixed-phase isotherms 
le and that T 1s 
identical with S, so that G always moves towards 5, 

Now 1s drawn to intersect at 
PW. ete. 
marked 


through 


GS, is drawn to imtersect 


G's gives the curve 


in Fig. 11. 


Joining up the 


operating line 


Determination of N, 


\ sloping line is drawn through the pomt | 
intersects the line M l. 


where the N scale 


mos 


Then the construction of section 2.5 ts carr d out, 


S replacing the T mentioned there, and the 


constructions of Fig. 9(a) or 9(b) being used as 
Details of the 


shown in Fig. 11. 


convenient, construction are 


\ vlance at 


is so nearly three 


Finally, the pont o, is obtained. 


the scale shows that N 


m 


vas side mass-transfer units of contact surtace 
are needed to accomplish the required cooling of 
the water. 


Comments on the calculation 


The value of N 


to within about 5 per cent for the calculation as 


, - can only be taken as accurate 
a larger number 
taken. 


is often good enough for 


described ; for higher accuracy 


of steps should have been However. 


5 per cent accuracy 
cooling-tower design. 

No numerical work was needed in the construc- 
tion. 


The make-up water rate, rw, Hiiw., 1S easily 


calculated from 


mG, 


NG, 


O-O161 
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the numbers being read off quite easily from the 
horizontal m scale. 


The ratio of inlet water to inlet air is given by 


Mw, 


m l Nip 
Pp 


Here the accuracy 1S poorer because of the 
cramping of the m scale on the left. 


Comparison with conventional calculation methods 
The data of the above example were chosen to 
accord with those of a problem worked out in 
detail by Treynat [9], who used the 
method [10], also with five temperatures intervals. 
TREYBAL’s results, with those of the present 
paper in brackets for comparison, were 


3-1 (2-0) 
IW 1-4 (1°37) 
(iw, Mwy) Me, O-O82 (0-026) 


Discrepancies therefore exist, Since the 
MerKEL method is in principle approximate, 
whereas the present one is so merely in practice, 
only a more detailed study would show which 
set of results is nearer the true answer ; however 
TREYBAL’s figure for the make-up water quantity 
is only an estimate. 

The times taken to determine the N.T.U. 


by each method are probably about the same. 


The present method shows to advantage however 
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when the moisture contents are sufliciently high 
to make the MerKken 
inaccurate. 


approximation too 


4, CONCLUSIONS 
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piece of transfer equipment may be carried out 
entirely graphically, apart from the final reading 
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Heat transfer in nozzles 


W. C. Racspate* and J. M. Smirut 


Purdue University, Lafayette, Indiana 


( Received 3 February 1959) 


Abstract—Local heat transfer rates and coeflicients were measured between steam and the 


wall over the converging and diverging sections of a nozzle having a diameter of 0-400in, at the 
throat. The coeflicients ranged from 43 to 290 B.t.u, (hr sq.ft) F) over a Mach number interval 


from 0-46 to 1-30. Maximum coefficients were observed at the nozzle throat. 


Correlations for well-developed turbulent flow in pipes could not be adapted to represent 
A simple flat-plate equation based upon a length Reynolds number pre- 


the observed results. 
dicted Stanton numbers about 10 per cent lower than the experimental values. Bartz’ method 


of modifying the flat plate concept to allow for gradients of bulk temperatures, velocities and VOI 
pressures along the nozzle axis agreed well with the experimental data and represented an im- ; 1 
provement over the flat plate equation. 

19590 


Résumé— Les vitesses d’echange thermique local et les coeflicients ont été mesurés entre la 


vapeur et la paroi sur les sections convergentes et divergentes d'un ajutage ayant un diamétre de 


0,400in. au col 
Les coefficients s‘étendent de 43 4 290 B.t.u (hr sq.ft°F) sur un nombre de Mach variant de 


0.46 & 1.80. Des coeflicients maxima sont observés au col. 


Des relations pour un écoulement turbulent bien établi, dans des conduits ne peuvent ¢tre 


adaptées a la représentation des résultats observés. Une équation linéaire basée sur plusieurs 


nombres de Reynolds prévoit des nombres de Stanton environ 10 fois plus petits que les valeurs 


expérimentales. 
La méthode de Bartz modifiant le concept linéaire pour tenir compte des gradients des 


températures globales, des vitesses et des pressions le long de axe de lajutage, conduit a un 


accord avec les résultats expérimentaux et constitue une amélioration de léquation linéaire. 


Zusammenfassung —Ortliche Warmeiibergangskoeflizienten wurden zwischen Dampf und der 


Wand in konvergenten und divergenten Teil einer Diise gemessen, deren engster Durchmesser 
ca. 10mm betrug. In einem Bereich der Machzahl von 0,46 bis 1,80 lagen die gemessenen 
Koeflizienten zwischen 244 und 1650 W m* grd mit einem Maximum im engsten Querschnitt. 


Die Beziehungen fiir voll entwickelte turbulente RohrstrOmung konnten die gemessenen 
Werte nicht wiedergeben. Die Gleichung fiir die ebene Wand mit einer auf die Lange bezogenen 


Reynolds-Zahl ergaben Stanton-Zahlen, die ungefiihr 10°, niedriger als die gemessenen Werte 


waren. Dagegen erheilt man eine gute Ubereinstimmung mit den gemessenen Werten, wenn man 
nach Bartz die Gleichung fiir die ebene Platte fiir Gradienten der Mischtemperatur, der Ge- 
schwindigkeiten und der Driicke lings der Diise verwendete; diese Methode stellt cine Verbesser- 


ung der Gleichung fiir die ebene Platte dar. 


INTEREST in jet propulsion has directed increasing nozzles is also unlike pipe flow in that the flow 
attention to the problem of predicting heat length is usually too short to establish fully- 
transfer rates in nozzles. Because of the large developed velocity profiles. It appears more 
gradients of pressure, velocity and temperature reasonable to describe the flow pattern in terms 
in the direction of flow, the conventional correla- of a constant-velocity core surrounded by a 
tions developed for flow in pipes or on flat plates boundary layer through which the velocity 
would not be expected to be applicable. Flow in decreases to zero at the wall. 


*Present address: Pratt & Whitney Aircraft Co., East Hartford, Conn, 
+Present address: Northwestern University, Evanston, Illinois. 
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SAUNDERS and 


CaLpER [1] measured heat 
transfer coeflicients in the diverging section of a 
nozzle by passing hot combustion gases (at about 
865 °C) through a nozzle made up of sections 
separated from each other by an air gap. It was 
concluded from the magnitude of the heat transfer 
coeflicients obtained that the boundary layer must 
be turbulent. Substantiating evidence was avail- 
able from velocity-profile measurements made at 
the exit of the nozzle. 

The theoretical development of heat transfer 
coeflicients for turbulent boundary layers is 
difficult and at present several assumptions have 
to be made to attain solutions of the momentum 
Bartz [2] has solved the 


problem by assuming, among other things, that : 


and energy equations. 


(1) A turbulent boundary layer exists through- 
out the nozzle. 


(2) The velocity and temperature difference 


(stagnation temperature — wall tempera- 
ture) variations in the boundary layer 


follow the 1/7 power law. 
(8) The friction in 


for a flat plate with the same boundary 
layer thickness, 


the nozzle is the same as 


(4) Reynolds analogy is valid for the boundary 
layer. 


SIBULKIN [3] solved the momentum and energy 
equations by that the 
layer was constant 
(incompressible boundary layer). Other theore- 
tical studies, for example that of KaLikuman [4] 
for the case of stream-wise pressure and velocity 
gradients, appear to be too complex to be useful 
in correlating heat transfer coeflicients. 


supposing, in addition, 


density in the boundary 


Sufficient experimental data are lacking for 
vareful evaluation of the theories of heat transfer 
in nozzles. 
not 


SAUNDERS and CaLpEr’s results do 
include the the throat 
the heat coellicients are at a 
maximum. Bopen [5] measured heat fluxes in a 
rocket motor fired with a nitric acid and aniline 
furfuryl system, but insuflicient data were taken 
to obtain heat transfer coefficients. Similarly, 
Zurcrow and Beighley [6] measured rates, but 


critical region of 


where transfer 


coefficients, in a rocket motor using a 
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nitric acid—jet-engine—fuel combustion mixture. 
GREENFIELD [7] employed a transcient method 
based upon the temperature-time history of a 
rocket motor copper. The results 
scattered considerably so that only an approxi- 


made from 


mate, empirical correlation could be attempted. 

The objective of the present work was to obtain 
reliable experimental data over the entire length 
of a nozzle, using a fluid whose physical properties 


were well established. Information was obtained 


for one nozzle design (throat diameter of 0-40 in.) 
at eleven positions between the entrance to the 
converging section and the exit of diverging 
region. 


APPARATUS 


AND METHODS 


Two stainless-steel, water-cooled nozzles of 


identical size and contour were used. Super- 


heated stream was employed as the heat transfer 
medium because its Prandtl number is close to 
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the 


recovery factor r has been found by a number of 


unity over a range of temperature. Since 


investigators [8-10] to be approximately equal to 


(Pr)'’*, this means that the adiabatic wall tempera- 


ture is equal to the more easily established stag- 
nation temperature. 

The the 
determined in both nozzles by measuring tempera- 
that 


heat transfer rate at wall, g,,. Was 
tures in the walls. It 
the 


wall. In 


re- 
the 


was 


Is necessary 
transfer rate normal to 
the 


approached by insulating various sections of the 


present heat 


nozzle sectional nozzle this 
nozzle from each other so that energy could flow 
in each section only in the radical direction. In 
the 


around the boundaries of a region inside the wall. 


solid nozzle temperatures were measured 
Then the temperatures within this region were 
obtained by solving the pertinent boundary value 
this 


profiles, and hence the radial heat transfer rate, 


problem. In manner radial temperature 
could be determined at any position along the 
nozzle length. 

The apparatus as a whole is schematically 
Fig. 1. Saturated 


was throttled to 5-20 p.s.i.g., 


pictured in steam passed 
through a strainer, 
and heated to the desired temperature level by 
an electric heater. The superheated steam then 


flowed through a straight, insulated, standard 


2 in. pipe to the nozzle entrance. In some runs, 
four 100-mesh screens were placed just upstream 
from the nozzle entrance to eliminate potential, 
thermal and momentum boundary layers at the 
beginning of the nozzle. 


Both 


steel (No. 416). Flanges and cooling jackets were 


nozzles were machined from stainless 
made from ordinary steel and copper wire was 
wrapped around the nozzles, inside the cooling 
jacket, to provide helical passageways for the 
and divergence 


cooling water. The convergence 


of the nozzles was slight, as shown in Fig. 2 and 
3, which give all the dimensions of the units. 

In the sectional nozzle, circular grooves were 
cut in the walls (Fig. 2) so that there were eleven 
sections of equal width and two end sections. Two 
ty in. thermocouple wells (one long and one short) 
were drilled in each of the eleven central sections. 
The grooves were evacuated in order to reduce 


the heat transfer from one section to another. 


and J. M. 

In the solid nozzle (Fig. 3), temperatures were 
measured at the same eleven longitudinal posi- 
tions. Two (one long and one short) ,& in. 
thermocouple wells were drilled in the nozzle wall 
at positions B, C, D, E, F, G, H, I and J corres- 
ponding to sections Nos. 2-10 in the first nozzle. 
Four the rmocouple wells were drilled (to four 
different depths) in the nozzle wall at positions 
A and K. 

The 
temperatures in 
O-O1 in, 


thermocouple used to measure 


the 


probe 


nozzle walls consisted of 
embedded in 
magnesium oxide inside a 0-062 in. 
O.D. stainless steel (No. 316) sheath. The short 


probes contained tron-constantan thermocouple 


two diameter wires, 


insulation, 


wires and the long probes copper-constantan. 


RANGE oF CONDITIONS 


Runs were made at four steam flow rates rang- 
110 to 200 hr. 


and pressures covered the range 516-617 “F and 


ing from Inlet temperatures 
5-20 p.s.Lg. The corresponding Reynolds numbers 
varied from 68,000 to 200,000, At each set of 
conditions four runs were made, two with screens 
at the nozzle inlet. and two without. 


CALCULATION OF Heat TRANSFER CORFFICIENTS 
Sectional nozzle 


From the temperature measurements, ft, and ¢,, 
at the two radial positions, the heat transfer rate 


for each section g, was obtained from the equation 


L (ty 


In ry 


t,) 


(1) 


where ¢ refers to the temperature measured at a 
An 


used to 


radius r from the centre line of the nozzle. 
the 
determine the wall temperature ¢ 
The 


(T,), was computed from 


expression of same form may be 


once g,, has 


stagnation 
the 


been computed. temperature 


following heat 


balance : 


(T inter 

In equation (2) the summation refers to the total 
the the 
axial position rc. The inlet stagnation temperature 


energy transferred from steam up to 
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was measured with a half-shielded thermocouple. Profiles, such as shown in Fig. 4, were extra- 
Since this is equal to the adiabatic wall tempera- polated to the diameter at the nozzle wall 
ture, T,,,,. for steam at the conditions employed, (designated by d,,). The slope and temperature 
the average heat transfer coefficient for the at the wall were then utilized in the following 


section is given by : equation to determine the heat transfer rate : 


dt 
} qu 8 Vu k oth 5 
(27 r,,L) (Tox te) (3) A dr (5) 
Solid nozzle As in the sectional nozzle, the stagnation tempera- 


- — ture at any location 2 was evaluated by an 
The temperature distribution in the walls of 2 : 

energy balance : 
the solid nozzle is described by the equation : 7 


dt »? A)dA 
t 1 2 t 0 (4) (Vu 
yr? 


( (T o)intet (6) 
for conditions of constant thermal conductivity. 
The boundary conditions, consisting of the Finally, the heat transfer coefficient was deter- 
measured temperatures indicated in Fig. 3, were mined from equation (3). 

such that equation (4) could not be solved 


ResuLts 
analytically. Relaxation methods were employed 


to obtain the complete temperature distribution, The heat transfer results are illustrated in 


using a Datatron digital computer (Electrodata Figs. 5 and 6*, where the local coefficient is 
Corporation). Illustrative results are shown in *Complete tabulation of the data and calculated results 
Fig. 4 for run no. 6. The letters identifying the re given in Ref. (12). 


radial profiles refer to the positions indicated in 


Colculoted volues , 


coefficient 


= 40} + + — 
§ 234567 89 WOH 


Section 


Fic. 5. Heat transfer coefficients for subsonic and 
supersonic flow of steam in a DeLaval nozzle. 
Sectional nozzle results. 


Predicted values : Experimental values 
Bartz with screens : 
flat plate, equation (11) © run 25 
flat plate, equation (10) @ run 36 
Average conditions : 
=e = 617 °F without screens : 


Fic. 4. Calculated radial temperature profiles at Pinter = 5 p-s-i-g run 29 
positions B through J for solid nozzle run 6. 1 112 1b hr @ run 40 
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Fic. 6 Heat transfer coeflicrents for and 
supersoni flow of steam in a DeLaval nozzk 
Solid nozzle results. 

Predicted values Experimental values 


Barrz with screens 
flat plate, equation (11) O run 
flat plate, equation (10) e run 16 


Average onditions 


lintet Without screens 
5 psig run 
Pinket 

ib he 2 


plotted versus the axial position in the nozzle. 
On each graph are shown the two sets of duplicate 


runs. one with entrance screens and one without. 


The coeflicients range from 3 to 2O0 
B.t.u ‘(hr ft? °F) over a Mach number variation 
from 0-46 to 18. The corresponding heat tlhux 


ranged from about 20,000 to 65,000 B.t.u (hr ft*). 

Fig. 5 indicates that the heat transfer coetlicient 
first decreases slightly and then increases in the 
converging section of the sectional nozzle. The 
throat. In 
the 


maximum value is attained at the 
the 


throat, the coetlicient rapidly decreases. Runs at 


supersonic region downstream from 
higher flow rates indicate larger coeflicients, but 
the variation with position is similar to that 
illustrated in Fig. 5. The average reproducibility 
of the data for the sectional nozzle as judged by 
the duplicate runs, using all the data, ts 3.5", 
for the runs 
screens. The results for the solid nozzle (Fig. 6) 


were not as good. Omitting the terminal positions 


with screens and without 


A and K, since these results were suspect because 


and 


of end effects, the reproducibility is 7-0°%, with 
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screens and 5-0°,, without. 
In the the heat 


coetlicients were approximately the same with 


sectional nozzle transfer 


and without entrance screens. In the solid nozzle, 


the results with screens are approximately 11°, 


lower than those without screens. While this 
seems to be somewhat larger than the repro- 
ducibility of the runs, it is difficult to explain 
how eliminating the thermal and momentum 


boundary lavers by inserting screens would 


the 
that the presence of the screens lead to erroneous 


reduce heat transfer coetlicient. It may be 
measurements of the entrance steam temperature, 
This comparison, along with the poorer reproduci- 
bility of the solid nozzle data, suggests that the 
measurements and computations for the sectional 
nozzle are more re liable. 

Comparision of the results for the two nozzles 
at the same flow rates (Figs. 5 and 6) shows that 
the the 


15 per cent higher than those in the sectional 


coetlicients im solid nozzle were about 
nozzle. An analysis of errors indicates that the 
data the the 
reliable. From the error analysis and reproduci- 
bility that the 


nozzle results are accurate within 10 per cent. The 


obtained for sectional nozzle are 


tests, it is believed sectional 


solid nozzle data showed more random variation 


and would be subject to larger constant errors. 


CORFFICIENTS 


Heat 
Flat plate methods 


PREDICTION OF PRANSFER 


\ simple concept of heat transfer in a nozzle 
that 


opened out, becoming a flat plate. 


is to imagine the circular wall surface is 
Because of 
the high Reynolds numbers in nozzle flow, the 


boundary layer is expected to be turbulent. 


Heat 


determined by 


transfer in this turbulent layer will be 
the the 


product of the fluctuating contributions of the 


time-average values of 
temperature and velocity. Suppose, in an approxi- 
mate fashion, the resultant effect is a velocity v’ 
Then 


the heat transfer rate per unit area through the 


perpendicular to the direction of flow. 
boundary layer ts 


Cy pi te) (7) 


: 

: 

a 

1959, 

‘3 

: 
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In terms of the heat transfer coeflicient this may 
be written : 
h=c, pv’ (8) 


St 


The effect of turbulence, as measured by v', would = 2 a 

be dependent upon the bulk velocity of the fluid © 
and the dimensionless ratio of the boundary 5 0-001+ 
layer thickness 5 and the distance x from the ®% 
leading edge of the plate. If v' is directly pro- 
portional to u and inversely proportional to 


Oo”, equation (8) becomes 


Length Reynolds No, (Re,)IO 


Fic. 8. Solid nozzle results (with screens). 


or Bartz method flat plate, equation (10) 


end value 


St a(; (9) 


Cp pu 


SaunpDERS and CaLper [1] used a modified 

The boundary layer thickness ratio is dependent form of equation (10) to correlate their data in 
upon the length Reynolds number Re,. Larzko the diverging section of a nozzle. The Reynolds 
[11] has developed an expression for this ratio number Re,. was based upon the distance 2’ 
» and then combined it with equation (9) to yield from the throat. The data of this study are 


the following equation : considerably lower than predicted by this means, 


as shown in Fig. 9. The SAUNDERS and CALDER 


St = 0.0292 (Re,) °? (10) 


expression may be written 


This expression along with the experimental 


St = 0.0285 (Re,.) °? (11) 


data are shown in Fig. 7 for the sectional nozzle, 


and in Fig. 8 for the solid nozzle. In evaluating 


the Reynolds and Stanton numbers, bulk values 


of the properties and velocity were used. Omitting é 


the data on Fig. 7 for the last section (No. 11) A 4 


in the nozzle, which are suspect because of end a= 


effects, equation (10) follows the trend of the 
observed results, but predicts values about 


10 per cent low. Though the results shown on 


Fig. 8 for the solid nozzle scatter more, the same 


conclusions apply. 


Diverging region results. 


end section value 


ay SAUNDERS and CALDER method, flat plate 
ad uation (11) 
experimental result, sectional nozzle 


© experimental result, solid nozzle (with screen) 


Nozzle equations 


The use of the flat plate approach (equation 10) 


Fic. 7. Sectional results (with screens). 


Bartz method flat plate, equation (10) neglects the possible effects of pressure and 
end section value velocity gradients, in the direction of flow, on the 
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heat transfer rate. Banrz [2] and SisucKin [3] 


have taken these gradients into account by 
integrating the equations of motion, again 


supposing that a turbulent boundary layer exists 
throughout the nozzle. 

Baxrz’ solution, based upon the assumptions 
given earlier, is in the form of two equations, one 
(equation 12) for the momentum boundary layer 
thickness, @, the 
temperature boundary layer thickness, A. The 
heat is then expressed in 


terms of ratios of 2 and A to the velocity boundary 


and one (equation 13) for 


transfer coetlicient 


laver thickness 4. 


54 
da s 


Me 2 /@) d 
| (4) 


d A 9) d [, AGc,(T, | 
dx dx | r 
0/8)(A/A 
fps i Pr® 46 
! 
Pr® 16 


14 
These three equations can be solved, numerically, 
for the heat transfer coetlicient and the boundary 
layer thickness at any point, from a knowledge 


of the 


the flow rate, wall temperature distribution and 


steam temperature entering the nozzle, 
the nozzle geometry. 

The Bartz equations do not indicate a unique 
Stanton 
number, Re, 


relationship between the number and 
length-Reynolds 


the results for different operating 


However. for 


each nozzle. 
conditions correspond to a single such relation- 
The 


dotted lines, indicate good agreement with the 


ship. results, shown on Figs. 7 and 8 as 
experimental data. 

SIBULKIN’s method {3}, based upon assuming 
the boundary layer IS incompressible and some- 
what different assumptions than Bartz made, 
was also employed to calculate heat transfer 


The results do not 


coellicients for the nozzles. 


and J. M. 


SMITH 


agree as well with the experimental data as the 


Barrz or flat-plate prediction methods. 


Application of pipe flow equations 


The Colburn equation for well-developed tur- 
bulent flow in pipes 


St = 0.023 (Re,)°? (Pry (15) 
(hte, 


was adapted to nozzles by using the point value 
of the nozzle diameter and the bulk velocity in 
The 


perties were evaluated at the bulk temperature 


formulating the Reynolds number. pro- 
of the steam. The Colburn equation is compared 
the the sectional 


screens) in The 


data for 
Fig. 10. 


with nozzle (with 


poor agreement in 


¢ 
= 
Ww) 
04 2 03 0 
-6 
Diameter Reynolds No, 


Fic. 10. 
Colburn equation, equation (15) 


Pipe flow equation (Colburn) correlation. 


r experimental data, sectional nozzle (with screens) 


denotes end section result 


comparison with Fig. 7 and 8 indicates that the 


boundary layer concept more appropriately 
describes heat transfer in nozzles. Other pipe-flow 
correlations [13, 14] were compared with the 
observed results, but no improvement was noted 


over the Colburn equation. 


Conc LUSION 


heat transfer coellicients in 


were found to decrease 


Experimental 
nozzles slightly at the 
entrance to the converging section and then 
the throat. 


In the supersonic, diverging section the coefticients 


increase to a maximum value at 
decreased sharply. 

The data agree surprisingly well with a simple 
equation based upon a turbulent boundary layer 


\ Ol 
{ 0-0288 | (12) 1959 
+ 
= 
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flowing along a flat plate. Better agreement is T, 


obtained with Barrz’ equations, which still 
presume a turbulent boundary layer but take be 
into account velocity and pressure gradients in 
the direction of flow. Prediction methods based u, 
upon modifying correlations for turbulent flow in v 
pipes are not as satisfactory. ” 
i 
NOMENCLATURF 
1 area ft? 
i. c¢ oss-sectional area of nozzle at throat ft? 
Cp Specific heat at constant pressure B.t.u/Ib 
d diameter of pipe or nozzle ft 
, G mass velocity Ib hr ft? 
h heat transfer coeflicient B.t.u hr ft? F 
’ hk. thermal conductivity of stainless steel 5° 
(416) (15-2 B.t.u/hr ft 
1] k = thermal conductivity of fluid in nozzle 
B.t.u hr ft °F 
L width of heat transfer aren in a single 
section of the nozzle ft 
V Mach number dimensionless A 
p static pressure p-S.1.g. 
Pr Prandtl number dimensionless 
q rate of heat transfer B.t.u hr 
’ radius of any point in the nozzle. r, and r, 
corresponds to radii at which temperatures were 
menusured ft 
The symbol r is also used to indicate the recovery \ 
factor, defined by the equation 
Tn 
p 
Rey Reynolds number dimensionless —— 
Pe 
Re, length Reynolds number dimensionless ~~? 
Subscripts 
T static temperature of fluid F 1 
adiabatic wall temperature F i 


REFERENCES 


Barrz D. R. Trans. Amer, Soc. Mech. Engrs. 1955 77 1235. 


4} Kanikuman L. E., N.A.C.A. Tech. Memo. No. 1229 1949. 

5) Bopen R. H. Trans. Amer. Soc. Mech. Engrs. 1951 73 385. 

Zocrow M. J. and C, M. Jet Propulsion 1952 22 323. 

8S. J, Aero. Sci. 1951 18 512. 

[8] Kaye J., Keenan J. H. and McApams W. H. Heat Transfer and 


{11} Larzxo Z. Angew. Math. Mech. 1921 1 268. 
{12} Racspate W. C. Ph.D. Thesis, Purdue University, August 1957. 


nozzles 


Saunpers O. A. and Catper P. H. Heat Transfer and Fluid Mechanics Institute 1951 p. 91. 


SipuLKin M, Heat Transfer and Fluid Mechanics Institute, Pape 


stagnation temperature of fluid F 
temperature within nozzle wall KF 
bulk or free-stream velocity ft, sec 
velocity in boundary layer, in the a 

direction ft sec 
velocity at nozzle throat ft se« 


velocity in y direction in boundary layer ft/sec 
flow rate of fluid Ibshr 
distance from entrance of nozzle measured in 
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defined by equation 
8 
5° | | dy 
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Abstract—An analytical solution 


tion and by the steady rotation of one or the 


capillaire ou rotatif, 


comment les équations peuvent tre 


INTRODUCTION 


IN RECENT years, a number of articles concerning 


the flow of non-Newtonian fluids in annular 


spaces have appeared. Votarovicn and 


[1], Sacurpanov [2] and van OLPHeEN [3] published 


approximate solutions for axial flow of Bingham 


plastic materials in annuli, and Mort and 
OroTakeE [4] published a more general solution 


Helical flow of an annular mass of visco-elastic fluid 
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of the equations of change is presented for the steady flow 


equations may be applied to other cases of interest. 


appliquées a d'autres cas intéressants 


gezeigt, wie die Gleichungen auf andere interessierende Falk 
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of an arbitrary viseo-elastic fluid through a concentric evlinder annulus. It is assumed that 
motion is imparted to the fluid by an impressed pressure gradient and or a gravitational accelera- 
other, or both, of the annulus cylinders. Under 
such conditions, the paths traced out by individual fluid particles will be circular helices 

The solution is presented in terms of definite integrals which contain an arbitrary function 
F(Y); this function may be called an “apparent viscosity.” It is shown how the function 
F(Y) may be determined from standard experiments with capillary, or rotational, viscometers. 


Suggestions are advanced for performing the required integrations, and it is shown how the 


Résumé —L auteur présente une solution analytique des équations de la variation, pour lécoule- 
ment permanent, dun fluide arbitraire visco-clastique, a travers un anneau entre deux cylindres 
concentriques. Il est supposé que le mouvement est communiqué au fluide par application 
dun gradient de pression ou par Paccélération de la pesanteur, ou par les deux actions simultanées, 
et par la rotation permanente de lun ou lautre, ou des deux evlindres constituant lanneau. 

Dans ces conditions, les trajectoires de chaque particule du fluide seront des hélices cireulaires. 

La solution est présentée en fonction d'intégrales détinies qui renferment une fonetion 
arbitraire F (Y). Cette fonction peut étre appelée “ viscosité apparente.” L’auteur montre 


comment la fonction F (¥) peut étre déterminée par des expériences standard avec viscosimétre 


Des suggestions sont présentées pour effectuer les intégrations nécessaires et il est montré 


Zusammenfassung —Es wird cine analytische Lésung fiir die Verformungsgleichungen einer 
stationiiren Stromung einer belicbigen viskoelastischen Fliissigkeit durch cinen konzentrischen 
zVlindrischen Ringspalt mitgeteilt. Die Bewegung kann der Fliissigkeit entwerder durch einen 
angelegren Druckgradienten und oder durch ein Beschleunigungsfeld und durch die stationiire 
Rotation des inneren, des Ausseren oder beider Zvlinder des Ringspalts aufgedriickt werden. 
Unter diesen Bedingungen sind die Stromlinien einzelner Fliissigkeitsteilchen Kre isspiralen 

Die Léasung wird in Form bestimmter Integrale mitgeteilt, die cine beliebige Funktion F (Y) 
enthalten; diese Funktion sei“ scheinbare Viskositit * genannt. Es wird gezeigt. wie die Funktion 
F(Y) aus Standardversuchen mit Kapillar- oder Rotationsviskosimetern bestimmt werden 


kann. Zur Ausfiihrung der verlangten Integrationen werden Vorschlige gemacht und es wird 


angewandt werden kOnnen. 


which unfortunately is in error. The exact 
solution for the flow of a Bingham plastic in an 
annulus was published by Larrp [5] and by 
FREDRICKSON and Birp [6]. Freprickson and 
Birp [6] also derived a solution for a so-called 
* power model ”’ fluid ; their results for both the 
Bingham plastic and the power model fluid were 


given in terms of dimensionless variables. The 
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have been extended by MELrose et al. [7] and by 
SAVINS [8]. 


FrepRIcCKSON [9] has derived a solution for 


the combined axial and tangential flow of an 
arbitrary, inelastic, non-Newtonian fluid in an 
annulus, His equations are actually a solution 
to a special case of the system of differential 
equations derived by Riviin [10] to describe 
the combined axial and tangential flow of an 
annular mass of visco-elastic fluid. It is the 
purpose of this paper to show how the equations 
of Freprickson may be generalized to give the 
solution of the problem posed by Riviin ; viz., 
for the helical flow of an annular mass of visco- 
elastic fluid. 


CONSTITUTIVE 
FLUIDs ; 


EQUATIONS FOR VISCO-ELASTIC 


Tureory or anp Ericksen 


Riviiw and Ericksen [11] have published a 
phenomenological theory of  stress-deformation 
relations for isotropic materials. On the basis of 
that theory, Riviiw [12] has shown that if one 
assumes that in a visco-elastic fluid which is 
isotropic in its state of rest, the stress components, 
7m, at a point in the fluid are expressible as 
polynomials in the gradients of velocity, accelera- 


tion, second acceleration, ..., (” — 1)* accelera- 
tion at that point and then the stress matrix, 
Nn ( \|7r45||) may be expressed as polynomial in 
n kinematic matrices D,, D,, ..., D,. The 


elements of the matrix D, ( D,||) are just 
twice the elements of the usual velocity-strain 


matrix : 


Di® = — + — (1) 
v2; 


and the matrices D, (= |)" ||) with r 2. 3. 
-» m are defined by the recursion formula 


A) 
dt i 1 
3 
+ 2 | D,,;" >) (2) 
mel oa; 


If the geometry of the flow situation is such 
that D,=0 for r > 2, then Rivirn [12] has 
shown that the relation between the stress 
matrix and the kinematic matrices is given by 


Helical flow of an annular mass of visco-clastic fluid 
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n pl D, Lo D, 43 D,? 
z, 
z;(D, D, + D, D,) + % (D,? D, 
D, D,*) 


x, (D, D2 ~ D,) 
(D,? D,? + D,? D,?) (3) 


for incompressible fluids. The coeflicients «, are 
not constants, but are expressible as polynomials 
in the ten* scalar invariants tr D,, tr D,?, tr D,’, 
tr D,, tr D,’, tr D,’, tr D, D,, tr D,? D,, tr D, D,?, 
tr D,? D,”. In the above expressions, it is to be 
understood that the usual rule [13] of matrix 
multipication is applied and that D," D/'D.. 
The trace (“ tr’) of a matrix is simply the sum 
of its diagonal elements. 

For inelastic fluids, the stress is a function 
only of the velocity-strain (i.e., up to an arbitrary 
hydrostatic pressure, 
D_ with r 


fluids, Rivirn’s equation becomes 


p), so that the matrices 
1 do not affect the stress. For such 


— pl (4) 
Kquation (4) is the equation derived by REINER 
[14, 15] and by Rivitn [16] and applied by Rrviin 
[16, 17] to explain the so-called “ Weissenberg 
effect [18]*. 
called the viscosity 
has been called [19] the “ 


The coefficient «, may thus be 
and the coefficient a 
normal stress coefficient.’ 
The derivations of Freprickson cited above 
utilized the constitutive equation (4) with the 
further hypothesis that the coeflicient a, could 
be set equal to zero. 

O_proyp [20] has formulated an alternate 
theory of visco-elasticity. The constitutive equa- 
tion of OLDROYD contains terms not only in the 
gradients of velocity, acceleration, ..., (n — 1)** 
acceleration, but also in the convective derivatives 
of the elements of Nl. 


it will be assumed that the rheological behaviour 


In this paper, however, 


of the fluids in question is given by equation (3). 


*For incompressible fluids, tr D, 0, so we need 
consider only nine scalar invariants. 
tOLDROYD [20] has sugvested that the Weissenberg 


effect may be due to fluid elasticity ; this view seems to 


be borne out by the experimental work of Roberts [21]. 
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E@uaTIONS FoR FLtow ANNULI 


Consider the steady, combined axial and 
tangential flow of a mass of visco—elastic fluid in 
a concentric cylinder annulus. Axial motion is 
imparted to the fluid by an impressed pressure 
gradient and/or a_ gravitational acceleration ; 
tangential motion is imparted by causing one or 
the other or of both the annulus cylinders to 
rotate with constant angular velocity. By 
hypothesis, we are dealing with a region far 
enough removed from the entrance (or exit) of 
the annulus so that the angular velocity, w, and 
the axial velocity, u, at any point of the region 
depend only on the radial distance of that point 


from the axis of the annulus. 


Let the z-axis of a rectangular Cartesian 
co-ordinate system lie along the axis of the 
annulus, and point in the direction of axial flow 
(Fig. 1). Then, if one chooses the direction of 
the a-axis so that y = 0 at the point of fluid 
considered, the kinematic matrices D, are given 


by [10]; 


| 0 re | | 
0 0}; D, 0 0 0); 
} u 0 0 | | 0 0 0) 
D.=@, 2 (5) 
where primes denote differentiation with respect 
to r. 
Let 


2 (r? + u’®) (6) 


Then the invariants of D, and D, may be written 
as {10}. 


trD,=0, trD,* 0, 

D,= ¥*, tD,* = Y*, 

tr D, D, = 0, trD? D, = 4 Y*, 

tr D, D = 0, and (7) 


Hence, it appears that all coeflicients a in equa- 
tion (3) are functions only of Y. But Y is a 
function only of r so equations (7) show that the 
x are functions only of r. 

Substitution of equations (5) into equation 


(3) yields the expressions for the stresses : 
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PATH OF FLUID 
PARTICLE 


to DIRECTION OF GRAVITY 


Fic. 1. Helical flow in an annulus 


p (22, %,) ¥ ¥? 
a, 
Tio rw’ [x +a, + a 
= [ x, + ag + a, 
Tes rw (3) 


However, by integration of the equations of 
motion, Riviin [10] showed that 
Br; ™13 Pr (9) 
where B and C are constants of integration and 


P (another constant) is defined by 


P p (10) 


In equation (10), p is the density of the thud 
and g, is the z-component of the gravitional 
acceleration g. (It is assumed that g, and g, 
may be neglected). The quantity P may be 
called the frictional pressure grachent. 

Define the dimensionless radial co-ordinate & by 


E=r/R (11) 
in which R is the radius of the outer annulus 


cylinder. Then combination of equations (8), (9) 


and (11) gives 


ae 
. 
Q, 
KR 
R \ x 
| 
‘ 
: VOL 
l] 
q 1959 
i 
= 


= = + a, } a, (12) 
PR? (& 
Rm = | 


du 
= Te + as 


where the constants 8 and A may be derived 
from B, C and R. Equations (12) and (13) are 
identical to equations (4.9) of Riviin’s paper 
[10]. 


Since the «, are functions of Y, we can set 


F(Y)=a, +a, Y +a, ¥? (14) 
so that (12) and (13) become 

dw 

dé (15) 


dé 


The separation of variables represented by 
equations (15) and (16) is permissible, since Y, 
and therefore F(Y), is a function of r(or &). 
The velocity distribution is then obtained by 
integration : 


PR? dt 
| 


F(Y) 


In equations (17) and (18), we have used the 
boundary conditions that 


w= $2 


at r= R(é 1) (19) 


where 2, is the angular velocity of the outer 
cylinder of the annulus. The boundary conditions 
at the inner cylinder of the annulus are 


ao => $2; 


at r K) (20) 
u=0 


where «R is the radius and 2, is the angular 
velocity of the inner annulus cylinder. The 


Helical flow of an annular mass of visco-clastic fluid 


a, (13) 


conditions (19) and (20) are the statements that 
the fluid does not “ slip si 


at the solid boundaries 
of the flow system. 


From equations (17), (18) and (20), there result 


1 
dv 
Q,— 2,= (17a) 
1 


Equations (17a) and (18a) are the determining 
equations for 8 and A. A method for performing 
the integrations called for by these equations 
will be presented below. 


The volumetric flow rate, Q, is 


OR 
u(r) rdr dO | (21) 


0 
or, because of (18): 
; 
Q | eae | (© 


. . > 


- (22) 
F(Y) 
Interchange of the order of integration and re- 
arrangement gives finally 


PR 


Equation (23) is the desired expression for the 
flow rate through the annulus. 

The values of the torques required to maintain 
the cylinders of the annulus in steady rotation 
may now be computed, The torque on the inner 
cylinder (per unit length of cylinder) is given by 
the expression : 


*iquations (17a), (18a) and (23) show that for a given 


fluid, the relations between 4Q rh, PR 2, « and 


9 22; are independent of the scale of the apparatus 


upon which those relations are determined. This would 
be an important consideration in any experimental 
programme designed to test the validity of the theory 


developed herein 
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dt 
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1 

7 — (72 — 

s= - (18) F(Y) 

= 
u=0 
= 
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or, because of (12), 


G R* B, 


(24) 


Since in a steady state or rotation, the same 
couples must act on every lamina of the fluid, 
it follows that equation (24) also gives the torque 
acting on the outer cylinder (per unit length). 
DETERMINATION OF CONSTANTS; INTEGRATION 
or E@uations (17a), (18a) and (23) 


the constants 8 and A must be 


trial-and-error 


In general, 


numerical inte- 


In order 


determined by 
gration of equations (17a) and (18a). 
to perform these integrations, it Is necessary to 
have a graph, table, or analytical expression 
describing the Fi(yY) with Y. It 


will be shown below how such a graph, table, ete. 


variation of 


may be prepared by analysis of viscometric data. 
For “* pseudoplastic ” fluids, the graph of F(Y) 
would have the general form of Fig. 4 in the 
article of and Roperts [22] |for 
(apparent) 

At any rate, suppose viscometric experiments 
have established the dependence of F Y) on 
Y. Then to perform the integrations, the following 
procedure 1s suggested. 


8 and A are assumed, 


\ alues of Ti 


Reasonable values of 
One then computes the and m4, 
from equations (12) and (13) for a series of values 
of & between & x and &€ 1. At a given 
radius, €, a value of F (Y) is assumed and rw’ 
The 
value of Y is then computed, and from the 
graph, table or equation for F(Y) vs. Y the 
value of F(Y) corresponding to the computed 
If the assumed value of 


and u’ are computed from equations (5). 


value of Y is found. 
F(Y) checks with the computed value of F(Y), 
one then proceeds to the next value of & in the 
series and F(Y) at that radius; 
otherwise, new values of F (Y) must be assumed 
and the calculations repeated until assumed and 
calculated values of F (Y) agree. 

Once the variation of F (Y) with & correspond- 
ing to the assumed values of 8 and A has been 
established by the process described above, the 


computes 


integrations called for by equations (17a) and 


*In that article, the abscissa is Y rate of shear 
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(18a) are performed by some numerical means. 
If the assumed values of 8 and A satisfy equations 
(17a) and (18a), and 
equations (17), (18) and (23) may be integrated 
(numerically) without difficulty. If the assumed 


values of 8 and A do not satisfy equations (17a) 


those values are correct 


and (18a), new values must be assumed and the 
entire process repeated. It is obvious that the 


amount of computational work necessary ts 


rather formidable, and the method is best handled 
by a digital computer. 

Finally, let it be noted that the values of 8 and 
A so computed will hold for a particular set of 
values of PR 2, « and Q,— 2, If any of 
these quantities are changed, then 8 and A will 


change also. 


DererMINATION oF F(Y) rrom Viscomerru 
Data 

The equations describing the flow of a visco- 
elastic fluid in a rotational (“ Couette’) type 
viscometer may be derived from the equations 
given in this paper by setting u = 0. Hence, 
for the Couette viscometer, the kinematic matrices 
become 

00 
D, rw 0 0); D, 0 00); D, = 0 (25) 
oO 0 OO r>?2 


The ten invariants of D, and D, are again given by 


equations (7), but the invariant Y assumes the 


simpler form 
(tia) 


For the stresses, we find 


p + (2a, + a) } 
(x,+ %) + a, 
Tie rw’ [a x, rw’ F(Y) 
Tos 0 (12a) 


Hence, the function F (Y) is simply 


F(Y)="" Couette. (26) 
rw 
The stress 7,, may be determined from the 
equations 


. 
Z 
. 
‘ 
tne 
: 
‘ 
oi 
- 


Helical flow of an annul 


= B (27) 


G, = 27 R? B (24a) 


and the measured torques on the cylinders*. The 
quantity rw’ is the usual rate of shear. and may 
be determined from experimental data by the 
standard means [2! 

Equation (26) shows that F(Y) is nothing 
more than the viscosity (apparent viscosity, really, 
since the usual Couette viscometer does not 
distinguish between elastic and inelastic lluids) 
as determined by the Couette viscometer. Hence 
a plot of viscosity, as given by equation (26) 
against twice the rate of shear squared | y) 
gives the plot of F(Y) vs. ¥ required in solving 
the problem of helical flow in an annulus. 

The equations describing the flow of a visco 
clastic fluid in a tube may be derived from the 
equations for helical flow in annuli by setting 
w=. Thus, the kinematic matrices may be 
written as 


OO 0 0 
“00 0 OO 


D 0. » 2 (28) 


for flow ina capillary tube. Because of equations 
(28), the ten invariants of D, and D, are again 


given by equations | 


Y is 


but here, the invariant 
2u* (6b) 


For the stresses. we find 


(x, + a) 
T 22 P a, Y, 
712 93 0 


[a +a, ¥?]—u' F(Y). (12b) 


Hence, the function F (Y) is given by 


F(y)=™ 


Capillary (29) 


*The experimental data must of course be corrected 


for end effects, non-uniformity of shear rate. ete. The 


hecessary corrections have been discussed by Toms [23 


i~? 


ir mass of visco-elastic fluid 


In the capillary tube viscometer, the rate of 
shear is not uniform (as, to a first approximation, 
it is in a Couette viscometer with a small gap) 
but is a maximum at the wall of the tube and 
zero at the centre of the tube. However. the 
shear rate at the wall may be determined from 
the Rapinowrrscu—Moonry equation [24, 25], 
and the corresponding shearing stress at the wall 
Is given by 


(30) 
wall 2 
Hence, the function F(Y) mia be determined 


from capillary viscometer data by plotting 


™3 (again, an apparent viscosity) vs. 
twice the rate of shear at the wall squared ( Y). 
The data must be corrected for entrance and 
exit losses, kinetic ene rev ¢ ffects. wall effects. 
etc. Methods of correction have been discussed 


by OLDRoyp [26] and by Freprickson [27]. 


Discussion Conc. SIONS 


It was possible to solve Rivitn’s equations 
for helical flow in annuli because of the relatively 
simple geometry of that type of flow. Thus, in 
helical flow. there are only two (u and rw’) 
non-vanishing, independent elements of the 
velocity-strain matrix, and these are related to 
the stress (7. and 7,,) by the same relation. 
There is the further consideration that for the 
case considered, the coeflicients x. in the con- 
stitutive equation (3) are functions of only one 
invariant quantity, Y. As shown above, these 
two circumstances make it possible to separate 
the variables in Rivurn’s equations and so to 
affect a solution. 

In more complicated flow situations (such as. 
say, flow through a diffuser, de velopment of flow 
In a pipe, ete.), the simplifications noted above 
would not result, and any attempt to arrive at 
an analytical solution would be faced with grave 
difliculties. 

In this respect, let it be noted that the type 
of tlow which prevails in the usual viscometric 
experiments (with capillary tubes or rotational 
viscometers) is not very general, in that the 
rheological coeflicients x, become functions of 
only one invariant quantity (Y). In a more 


i 
4 
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general flow situation, the coeflicients x, would 


appear as functions of other invariants. Hence, 


the usual data obtained from capillary tube or 


viscometers may be inadequate to 


the 
Newtonian fluids in fairly general flow situations, 


rotational 


predict rheological behaviour of  non- 

The method of solution adopted herein may 
also be used to solve the problem of combined 
axial and tangential flow of Bingham plastic 
materials in concentric cylinder annul. According 
to OLDRoYD [28], the Bingham plastic ts deseribed 


by the constitutive equation 
n 
D, < (31) 
with 


(32) 


In the 


the (constant) plastic viscosity, and the mvariant 


equations, T, IS the ld stress, fig IS 


quantity 7? is detined by 


(33) 


For helical tlow of a Bingham plasty material 
im an annulus, one can easily solve the equations 


of motion and so obtam expressions for the 


stresses (equation %) and hence for the vield 


condition (7* 


. this will be a polynomial 


rin rt 


of sixth degree in the reduced radius 


solution of this polynomial for the two 
lving between « and 1 will give the boundaries 


of the The 


equations (31) may then be combined with the 


‘plug tlow region. constitutive 
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expression for the stresses and integrated to give 
the velocity distribution. 


Not ATION 


constants of integration 

kinematic matrix (elements 

torques acting on unit lengths of inner and outer 
evlinders of annulus 

acceleration due to gravity (components Z,, By: 

function of Y (“ apparent viscosity ") — defined 
by equation (14) 

identity matrix (clements 
frictional pressure gradient — defined by 
equation (10) 

hydrostatic pressure 

volumetric flow rate 

radius of outer annulus cylinder 

radial co-ordinate 

axial velocity 

i-component of velocity 

cartesian co-ordinates 

invariant quantity — defined by equation (6) 
th rheological coeflicient in equation (3) 
constants of integration 

Kronecker delta ( 


dummy variable of integration 


lifi=j; = ts )) 


‘ viscosity of Bingham plastic material 
detined by equation (32) 

ratio of radius of inner annulus cylinder to 
radius of outer annulus cylinder 

plastic viscosity 

stress matrix (elements mj) 

dimensionless radial co-ordinate 

density of fluid 

detined by equation (35) 

vield stress of Bingham plastic material 
angular velocity of inner annulus cylinder 
angular velocity of outer annulus cylinder 


angular velocity of fluid at 2, y, 
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This paper deals with unsteady state heat transfer im fixed beds in the most general 


Abstract 
The gas and 


ease of non-uniform initial bed temperature and varving gas inlet temperature 


described in terms of the solution to the more simple problem determined 


bed temperatures are 
initial bed temperature and constant gas inlet temperature, which solution has 


by uniform 
already been evaluated in several different ways 

|. leaving open the choice of the most practi al form of the 
the ‘ 


The treatment is kept genera 


solution to the simpler problem. In applications this choice is made, taking into account 


nature of the specific problem at hand. It will depend on the ranges of the independent variables 


Most fixed bed applications will be in a range where the 


error function solutions 


encountered 


aceording to KLINKENBERG are adequate 


non-stathonnaire 


Résumé —La communication traite de la transmission de chaleur a état 


dans des lits fixes Elle déerit le cas le 


plus general celui dune temperature initt ile non-uniforme 


variable du gaz a Tentree Les tk mperatures du gaz et du 


c'est a dire de celui a temperature 


dans le solide et dune temperature 
sont dérivées de la solution du probléme le plus simple 
et & température constante du gaz a lentree. Cette 


solution a 


initiale uniforme dans le solide 


déja été évaluée de plusieurs fagons lifferentes 


une liberté de choix dans le mode de solution 


Le traitement est tout-a fait veneéral, laissant 
ipplications pratiques, ce choix se fait en tenant compte 


du probleme le plus siniple Pour les 
li déependra des regions des variables indépendantes qui se ntent 


de leur nature speécitique 


‘tions de lit fixe les solutions a fonction de probabilite se lon KLINKENBERG 


Pour la plupart des applic: 


seront adequates 


Zusammenfassung Das Referat betasst su h mit der nicht-stationiren Wirmetibertragung 


vleichmassigen Anfangsbettemperatur und 


in Festbetten im allgemeinsten Falle der ni ht 
Bettemperaturen werden dargestellt an Hand 
Anfangsbettem- 


variabelen Gaseintrittemperatur Gas- und 
der Lésung des einfacheren Problems das bestimmt wird durch gleichmiissige 
Jsung bereits in mehrerlei verschiedener 


peratur und konstante Gaseintrittemperatur welche I 


Weise entwickelt worden ist 


Die Behandlung wird allgemem 


vehalten, wobei die Wahl der praktis« hsten Form der Losung 


Bei den Anwendungen wird diese Wahl ge troffen 


des einfacheren Problems offengelassen wird 


unter Beriicksichtigung der besonderen Art des jeweiligen Problems. Sie rik htet sich nach dem 


Die meisten Anwendungen von Festbetten 


Bereich der vorkommenden unabhangigen Variabelen 
Bereich liegen, in welchem die Fehlerfunktionslésungen nach KLINKENBERG 


werden im emem 


ausreichen. 


simple initial and boundary conditions as ~ the 


IN1 RODUCTION 
clementary problem.” A survey of the various 


Many authors working in the field of the heating 


or cooling of a packed bed by means of a flow 


forms of its solution, together with an examina- 


of gas have considered a bed at initially uniform — tion of their convenience of application and, for 


gas approximate formulations, of their accuracy as 


temperature being heated or cooled by a gv 
stream of constant inlet temperature*. We shall well, has been given by one of the present authors 


refer to the heat transfer problem with these [1]. 


*The usual assumptions include uniform gas and s id temperature over any bed cross-section, infinite thermal 


and absence of longitudinal mixing and conduction. 


diffusivity in the solid 
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A more complex case is that of a bed with an 
arbitrary initial temperature distribution, heated 
by a gas the inlet temperature of which is an 
arbitrary function of time. This problem has 
been treated by AMUNDson [2] and, more recently, 
by Reii.y [3]. Amunpson used Bessel functions ; 
Reit_y was the first to apply Fourier integrals 
or Fourier series. SAULBERG [4] recently published 


a solution to the 


mathematically identical 


problem of the steady state behaviour of a 
cross-flow recuperator in which the temperatures 
of the entering gas streams vary along the widths 
of their entrances. He used the same approach 
as developed on a more general basis in the 
present Bessel 


function solution to the elementary problem. 


article, and he employed the 

Nussevt [5] and Hausen [6] have solved the 
problem of the regenerator with an arbitrary 
initial temperature distribution, cooled or heated 
by a gas of constant inlet temperature. Nussevr 
Bessel 


Havusen’s Wirmepolmethode,” a heat source 


gave the correct function solution. 
method, resembles the approach of the present 
paper. Because of the symmetry of the governing 
differential equations in terms of the two depen- 
dent and the two independent variables, the 
treatment by HAvUsEN can 


easily be extended to include the case of variable 


Ni SSELT and by 


inlet gas temperature as well. 

We should like to point out that the analytical 
solution to the above general problem can easily 
be written in terms of any solution to the ele- 
mentary problem, without it being necessary to 
specify the form of this solution. It has already 
been shown in [1] that the solution to the ele- 
mentary problem can be expressed in a variety 
of formulations. Our treatment has therefore 
been kept general, in order to render possible the 
numerical solution of the general problem by way 
of any mathematical approach available for the 
elementary one, i.e. by any method yielding the 
response to a step-function input, The Bessel 
function solutions and our previous solutions 
employing error functions with corrected argu- 
REILLY’s 
method [3] is not suitable for this purpose, 


ment, for instance, may be employed. 


because the use of Fourier synthesis makes the 


evaluation of the step-function response super- 
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fluous. Although that response can easily be 
derived from his formulae, the use of it in the 
sense of the present article would be unnecessarily 
circuitous, 

In the following sections the derivation is 
given of two pairs of solutions to the general 
problem, expressed in terms of a function F. 
This function represents the total exchange of 
heat in the elementary problem as a function of 
bed height and time. It need not be known 


explicitly for this derivation, but its properties, 


as they are required, are given. Finally various 
specific forms of F which. by substitution, 


enable direct numerical solution of a variety of 
practical problems are discussed. 
Equations: Tue ELewentary 
Prosi EM 

The mathematical treatment of transient heat 
transfer in packed beds as proposed in this paper 
is subject to a number of simplifying assumptions, 
a few of which were mentioned in the footnote 
on p. 260, Many previous articles, such as that by 
Reintty [3], give a full account of them. 


The basic differential equations are 


T, T, 


\Z (1) 
(2) 
in which 
T, gas temperature 


solid temperature 


y dimensionless parameter proportional 
to length 


Z dimensionless parameter proportional 
to time*, 


In the definition of Z the time the gas takes 


to pass through the bed has been ignored, this 


*If l heat transfer coeflicient. a surface area of 


packing per unit volume, d = density, ¢ specific heat, 
t mean linear fluid velocity in the bed, F 


volume, 1 bed 


fraction by 
height, / time, index f fluid, 


index s solid : 
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time being very small in respect of the actual 
heating times. 

Equations (1) and (2) were solved by most 
previous authors for the conditions : 


Y=0-f, 1(Z > 0) 
(3a) 
Z=0-+T,=0(Y > 0) 
These refer to “ the elementary problem” of a 


bed initially at zero temperature, heated by a gas 
continuously entering at a temperature of unity. 


Heat 
Bep 


Tue Amount oF EXCHANGED IN 

THE 

The heat balance requires the heat lost by the 

gas to be equal to the heat gained by the solid. 

When written in dimensionless form this becomes : 
z 

(T,)y | (T,), (4) 


0 


F(Y,Z)=2Z 


The function F of Y and Z accordingly is the 
dimensionless representation of the amount of 
the amount of heat exchanged in the bed. It has 
been employed as such by San perc 

It is important to note that F is a symmetrical 
function, viz. 


F(Y,Z) = F(Z, Y) (5) 


This can be demonstrated as follows. Equations 
(1), (2) and (3) remain unchanged if the following 


pairs of variables are interchanged : 
Y and Z 
T, and 1 T, 
If the same substitutions are made in (4) we 


obtain : 


T, dl 


3rd member are identical 
(4), the 
relation (5), 


Since the 2nd and the 


with the 3rd and the 2nd member in 


symmetry of the function F, i.e. 


has been proved. 
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According to (4) the derivatives of F are closely 
related to 7, and T,: 


. (6a, b) 
)Z 
The relation between (6a, b) and (1) is obvious. 
As the conditions (3a) represent a unit step 
disturbance in gas inlet temperature, we may 
accordingly state : 
MF 
response of gas temperature to a 
unit step disturbance in inlet gas 
temperature. 
MF 
. response of solid temperature to a 
unit step disturbance in inlet gas 


temperature. 


In order to examine the effect of a unit step 
disturbance in the initial solid temperature we 
must consider equations (1) and (2) with con- 


ditions 


—0(Z> 0) 
(3b) 
Z-0>T, 1(¥ > 0) 


These conditions can be derived from (3a) by in- 
terchanging Y with Z and 7, with T,. 
transposition leaves the basic differential equa- 


The same 
tions unaltered. Hence the solution to (1) and 
(2) with conditions (3b) can be derived from the 
solution to the previous case by making the above 
substitutions. F being invariant with respect to 
and Z, this yields the 


following responses to unit step disturbances in 


the transposition of Y 


inlet solid temperature : 


MF 
1 . = response of solid temperature to a 
\ 
: unit step disturbance in initial solid 
temperature 
- response of gas temperature to a 
| 


unit step disturbance in initial solid 
temperature. 


All equations being linear, differentiation of a 
response to a step disturbance in gas or solid 
with respect to Z and Y, respectively, gives the 


response to a peak disturbance. Consequently 
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Unsteady 
the following peak response functions 
tabulated : 

to peak 
disturbance in: 


function : response of: 


gas temp. gas 
F solid temp. gas 
YF solid temp. solid 
gas temp. solid 


For most of the following discussion the type 
of the function F need not be specified, knowledge 


of the following of its properties being suflicient : 


(by substitution of Ga, b in 2 
secondly : 0; | 0 (8a, b) 
M4 0.Z Y.0 
(by substitution of 6a, b in 3a) 
thirdly : 1—e-’; 
IZ Y,0 
| . = ] e~* (9a, b) 
0,Z 


Equation (9a) is obtained by considering (7) 


along the boundary Z — 0 where, in view of 
(8b), (7) becomes a simple differential equation in 


dF DZ. 


manner. 


Equation (9b) is obtained in a similar 


Equation (9a) represents the well-known fact 
that the first gas which enters the bed cools down 
according to an exponential law because, during 
its passage, it constantly meets fresh solid at 
temperature 7, — 0. Similarly, the solid material 
at the entrance to the bed is heated by the gas 
according to an exponential law, since the gas at 
that place is constantly at the 
T, 1. This is expressed by (9b). 

Equations (7) to (9a, b) and their derivatives 


temperature 


enable considerable simplifications of the final 
formulae to be derived in the following section. 


SOLUTION OF GENERAL PROBLEM 


The (8a), 
representative of the elementary problem, having 


initial and boundary conditions 


served to determine the properties of the function 
F, will now be replaced by conditions (10). These 


determine the more general heat transfer problem 


state heat transfer in stationary 


can be 


2638 


packed beds 


subject to arbitrary gas inlet temperature and 


initial solid temperature. These arbitrary 

functions are designated by @. The new initial 
and boundary conditions are : 
Y=0+T, =80,(Z) 

(10) 
Z=0>T,=0,(Y) 


The problem of integrating equations (1) and 
(2) with conditions (10) can readily be solved 
without any further knowledge of the function 
F, by 


(HAUSEN’s method of heat poles 


superposition of responses to peaks 
The gas, 
for instance, passing the place Y at the time Z 
has the temperature which it would have had if 
no other disturbances had entered the bed before 

namely e* @, (Z) the 


grated contributions of all gas temperature peaks 


increased by inte- 
entering between the times 0 and Z and of all 
solid temperature peaks originally present between 
the places 0 and Y. The solid temperature is 


made similar 


up in a way. This vields the 
formulae : 
Zz 
0 
@, dy (lla) 
T,(¥.Z) 70,(Y) 
Zz 
| (sya) (11b) 


The reader wishing to check the correctness of 
these solutions by substitution in (1), (2) and 
(10) will find that he needs only the properties of 
the function F (7), 
(8a, b) and (9a, b) (including some of the deriva- 


expressed by equations 
tives of these equations) and not its specific form. 

In all previous work step functions rather than 
peak functions have been studied. Formulae 
(lla, b) may be adjusted accordingly by partial 
This has the effect of 


integration. replacing 
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integration over a series of successive 


magnitude @ by integration over a senes ol 


successive jumps equal to the first derivative @. 
the 


If use is made of equations (7) to (9a, b) 


resulting formulae become 


T,(¥. Z) = 0, (2) 


8, (0) 


,(0) — A,(0) } 


Psy), 


7 


Seprcipic Forms or F DiscUSSION 


For practi al appli ations the function F may 


be written in. for instance, the following specith 


forms : 


(15) 
JA 

It can easily be 
obey conditions (7) 
Y and Z is obvious. 

Substitution of (13) in (Ga, b) vields the 
mad 7, 
[10] (see 
Substitution of (13) in (lla, b) 


the ssel thon solutions whi hy were 


these forms 


proved that 
to (9a, b). The svmmetry in 
well 


known elementary solutions for 7, first 
wiven by 

survey |1}). vields 
che rived 


His 


solutions may be reduced to our simpler form by 


by Amunpsown ([2], his formulae 8 and 19). 


making use of the identity 


ye aks of 


HARMENS 


1, (2y/Z) dy (16) 


Nusseir [5 


though only for on 


(lla, b), 


boundary condition, 


obtained our solutions 


variabl 
his 


as inlet temperature being assumed to be 


constant (@, 0). 
Substitution of the 


first 


(14) im (Gab) vields 


elementary double power series solutions 


published by Sarria [7 
Formula (15) in (Ga, b) generates the solutions 


nm exponentials and double power series viven 
previously by one of the present authors | 1] 

The formulae (lla, b) 
(I2a, b) do not function F 
but This 


approximat 


two sets of final and 


contain the as such 


only its derivatives allows of the 


utilization of som solutions to the 
elementary problem which, by virtue of (6a, b). 


vield MF and dZ. 


easily bye 


although these forms 


cannot integrated and do not exactly 


obev the conditions (7) to (a, 


For 
the derivatives of F 


and Z 


approximated by 


all except very low values of y 


may be 


jerf{y Z 8/7 (17a) 


or, with only a slight loss of accuracy, by 


ly (Z 


as introduced by KuNKENBERG [1,! 
that 17a. b) 
longer of engineering accuracy for y 


that 


been shown l solutions 


and or 


below 2, but they rapidly become more 


variables 
Z of an S-shaped break 


the error was determined at about 


accurate for higher values of these 
Im the 
through curve 


1 Vy 


As a more accurate approximation, especially 


midpoint 


for low values of the variables, ONsacen’s 


| 
: 
’ 
WAR 
) 
‘ 
1959 
F= | 1, (2 y/ nb) dy d 
| 
berf | y(Z yi) (18b) a 
It has 
E 


solutions [1] without the added series may 


serve. In the present notation these are: 


berf Z \ 
VAL 
gist’ (19a) 
MF 
ert Z 
Z r \\ \ ) 
gia *1,(2/¥Z)  (19b) 


\ comparison of the preceding formulae leaves 
no doubt that (17a. b) and (18a. b) are the 
simplest ones for the calculation of the integrands 
for graphical or numerical integration of (lla, b) 
or (l2a, b). They can be used for all except 
very low values of Y and or Z which, however, 
in most cases imposes no practical limitation. 
The method can be applied in a much shorter time 
than Fourier synthesis, For problems 
with low values of the variables (such as with 
cross-flow recuperators ) the derivatives of 
functions (13), (14) or (15) or Onsacer’s solutions 
(19a, b) can be used for the evaluation of the 
integrands. 

It may be recalled at this point that elaborate 
tabulations of JF )¥Y have been prepared by 
BRINKLEY [9] for values of Y and Z between 0 
and 500. The other derivative, )F )Z, can also 
be read from the tables since F is a symmetrical 
function. Accurate as these tables are (six figures 
behind the decimal point), their usefulness is 
greatly impaired by the fact that the increments 
in the values of the independent variables are 
too large. This calls for non-linear interpolation 


in intervals of four or five figures, which is both 


cumbersome and inaccurate. 


Discussion or 


Returning once more to Rem.y’s paper [3] it 
must, unfortunately, be stated that his numerical 
example was illchosen to demonstrate the posst- 
bilities of his method. The initial disturbance in 
the solid temperature (his Fig. 1) can easily be 
described by superposition of two step-function 
disturbances. The problem may then be solved 
with the aid of the elementary solutions (6a, b) and 


the error functions with corrected argument 
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(18a, b), without any numerical integration being 
required, 

In example the solid temperature 
at any place in the bed can be described as 
1200 F plus the influence of a step disturbance of 

150 F plus the influence of a step disturbance 
of 100 F which has travelled over a slightly 
shorter distance : 


T, = 1200 150 ( -} + 100 | - | F (20) 
), is the dimensionless distance between the 
measuring point and the beginning of the bed; 
Y, is the dimensionless notation of a distance 
which is 0-5 ft shorter. Formula (18a) for large 
values of Y and Z can be written in a slightly 
simplified form : 
b+ — 4) —7/Y)] (21) 
In order to check Retiy’s curve for 7. after 


25min (this time being characterized by 


Z 97-22) a simple tabulation of )F DY is 
prepared according to (21). For the sake of 


convenience the points are taken 0-5 ft apart. 
Their distance from the bed entrance is designated 
by in ft. The following values of ean 


easily be derived : 


} 18-6 1-0000 
2 64-8 0.9942 
2-5 S1-0 0-83811 
3 97-2 0.4864 
3-5 113-4 0.1256 
rT 129.6 0-0142 
1-5 145-8 0-0007 


By proper substitution of )F.)¥ in (20) the 


following results are obtained : 


0-9942) + (100 1-0000) 


1150-9 F 


T, = 1200 — (150 


O-SS11) + (100 0-9942) 
1167-3 °F ete. 


T, 1200 — (150 


aft T, 1215-1°F 
3-5 1229-8 
1210-4 
4-5 1201-3 
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These figures are considered accurate to within the last decimal given. They check very well with 
£ 


the curve in Reiiiy’s Fig. 1. 
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The viscosity of liquids as a function of temperature 
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Abstract—In the past various formulae have been put forward as descriptions of the way in 
which the viscosity of liquids changes with temperature. Those proposed by Guzman and 
ANDRADI 

log 
and by SoupERs 

low T BR 


do not appear to possess any general validity. Of the formulae that have been published since 


then, that of Cornetissen and WarerMan 
log B 


is a good description of viscosity as function of temperature. The derivation of this formula 
from that of SouprErs is discussed. The one is evolved from the other by way of viscosity values 
for alkenes and other hydrocarbons, the ConneLissen WarerMan formula also being applied to 


ulass. 


Résumé — Différentes formules ont déja éte établies pour représenter les variations de la viscosite 


des liquides avec la température 


Celles proposées par GusMAN et ANDRADI 


et par Soupers 


log 


ne paraissent pas posséder une validité générale a formule publiée ultérieurement par 
CORNELISSEN et WATERMAN 


log B 


donne une bonne représentation de la variation de la viscosité avec la température. Létablissement 
de cette formule a partir de celle de Souprrs est discutée. Loune est dégagée de Pautre au moyen 
des valeurs de la viscosité des alkénes et d'autres hydrocarbures : la formule de CORNELISSEN 


WATERMAN s‘applique aussi au verre. 


Zusammenfassung In den letzten Jahren sind verschiedene Formeln fiir die Anderung der 
Viskositaét von Fliissigkeiten und der Temperatur vorgeschlagen worden. Nach Vorschligen 
von GuZMAN und ANDRADE gilt 

log » 
und nach Soupers gilt: 


RBeide Gleichungen scheinen keine allgemeine Giltigkeit zu besitzen. Von den Formetn, die seitdem 


veréffentlicht wurden, stellt die von CorNnetissen und WarrerMan 
log » 177. RB 


gute Beschreibung des funktionellen Zusammenhanges von Viskositét und Temperatur dar. 

Die Abswichung dieser P'ormel von der von Soupers wird diskutiert, Die eine wird aus der 
anderen aus Viskosititswerten fiir Alkene und andere Kohlenwasserstoffe entwickelt ; daneben 
wird die Cornnetissen WarerMan-Formel auf Glas angewandt 
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INTRODUCTION 


Various writers have proposed relationships 


between the viscosity of liquids and their tem- 
In 1930 AnpraApE published one that 


peratures. 
had been proposed by GuzMAN as far back as 


The appropriate formula ts 


log » AT+B8B 


1913 [1]. 
(1) 


viscosity, T the 


and Af 


the dynamic 


absolute temperature in degrees Kelvin 


in which » is 


and B constants. 


formula results in large 
discrepancies this 


revised the relationship in 1934, publishing the 


In many cases this 


and for reason ANDRADI 


following formula : 
log (» v' (2) 
Where 


it is a matter of displaying viscosity changes with 


in which v is the specific volume in @ mil. 


temperature in graphical form, formula (2) ts 
difficult 
1937 


to work with. Soupers [2] simplified 


it in into 


J. M. Ste 


vets and H. I. WarrerMan 


log v A/T +B (3) 


where v (kinematic viscosity) 7» p and p is the 


density in g ml. 

This formula was stated by the author to be 
just as satisfactory as (2). 

Since 1937 many other viscosity formulae have 
been proposed including that of CoRNELISSEN 


and WarerMan [3], which is 


+ 8 (4) 


in which A. B and x are constants. 


By quoting extensive numerical data, 
CORNELISSEN and WarerMAN demonstrated that 
their formula could be applied satisfactorily to 
very different 


and mixtures thereof, saturated and unsaturated 


systems, such as hydrocarbons 


mineral oils, and solutions of sugars and other 
substances. 
We 


between 


the relation 


SOUDERS In 


to discuss 
that of 


In order to give an impression 


now vo on 


shall 

this 
greater detail. 
of the magnitude of the errors arising when the 


formula and 


4 
| 
| 
| 
tcrnene 
Pa | 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
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The viscosity of liquids as a function of temperature 


CuzMAN-ANDRADE and Soupenrs formula are em- 


undecene 


ployed, the viscosities of ethylene, 1 
plotted 
temperature [4] in Figs. 1, 2, 3 and 4. 


and l-eicosene have been against 


In Figs. 1 and 2 are represented not only the 


values of the logarithm of the viscosity but also 


the limits of the experimental accuracy by two 
horizontal lines. 

In Figs. 3 and 4, however, the experiments are 
so accurate that the points sufficiently represent 
the experimental values. 
that the 


The graphs in Figs. 1 and 2 show 


thene if 
| 
| 
VOL, 
lL] 
Fic. 2 
0-60} 0-70} 
ene -Eicosene 
050} 60} 
030} 0-40} “a 
g 020) 030} 
O10} O20} 
Undecene 
-0 30} -0-20} 
| 
; 26 27 28 29 30 31 32 33 34 26 27 28 29 3C 3 32 33 34 
1000 / T 1000 /T 
i Fic. 3. Fic. 4. 
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Table 1. 


WATERMAN 


undecene 


©) E. 
(keal 1) 


(keal mol) 


1 crcoscene 


| 
(keal mol) (keal mol) 


P20 


7s 


GuzMan-ANDRADE and Souprers formulae have 


a reasonable degree of validity for ethvlene, 
if we take the expe rimental accuracy in considera 
For 1 and l 


pancies (see Figs. 3 and 4) are so great 


discre 


that it 


tion. undecene eiosene the 


is qpuite impossible to apply the formulae (1) and 


(3) on these Therefore straight dotted 


lines are drawn which show the deviations. 


Their failure is connected with the fact that 
d(log») T) and dilogyv) dil T are not 
constant but are dependent on temperature. 


In various places in the literature [5] the term 


activation energy is used to indicate the value of 
this quantity multiplied by 2-30 R (# is the gas 
1-O87 cal Hence 

forward we shall be using EF. and E, to stand for 
2-30 Rd(logy) and 2-30 Rd(logv) dT) 


respect 1\ ely 


constant. whose value Is mol 


Table 1 contains figures showing how Ey and 
and Ev depend on temperature. 
The Table reveals that FE, and EF 


decrease being all the 


decrease as 
temperature increases, the 
more rapid according as the molecular weight ts 
If the dependence on temperature ts 
T*'. in which ¢ 


formula of CORNELISSEN 


greater. 
expressed in the form E, = ¢ 


and w are constants, the 


and WarrerMan can be derived from it by inte- 
gration. 
The relationship log» = A 7” B can be 


deduced from the above in an analogous manner. 


The value of « for hydrocarbons does not seem to 
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change verv much if the cde pendence on tempera- 


of log » is studied instead of that of log v. 


ture 


However, this is not generally true of all sub- 
stances, 
It depends on the way in which the specitic 


affected by 


temperature. For ¢ xamiple, in the case of mercury, 


weights of various substances are 


[3] it is found that x is about 1-5 when the log: 
formula is employed and about unity when the 
log » one Is ¢ mploved. The authors prefer to use 
the log » formula for liquid hydrocarbons because 
this the 


measurement, 


obtained from 


value is one actually 


Even so, the log » formula is sometimes useful 
when data given in the literature are being worked 
the » of various types of glass is 
different 
their specific gravity 1s not. 

1 he 


is applicable to simple hydrocarbons whether of 


on. For xample 


often known for temperatures, while 


formula of CORNELISSEN and WATERMAN 


weight and also to 


weight; it is 


large or small molecular 


mixtures of hydrocarbons. value of a 


increases with molecular about 


unity for ethylene and increases as the lengths 
of the chain does. «x for 1l-eicosene ts about 2. 


DereNDENcE OF Ey ann Ev on Cuain Lenoru 

Table 2 displays EF, and EF, values for various 
alkenes as determined at two temperatures, 
and 


the length of the chain. The same is found to be 


Table 2 shows that increase with 


VOL 


}OSGs 


' 
q 
2 2-22 
= 
q 
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Table 2 


(keal mol) 


heptene 
octene 
nonene 


decene 


undecene 3-00 
dodecene 3-22 
tridecene B46 
tetradecene 3-78 
pentadecene 


hexadecene 
| heptadecene 
1 octadecene 
nonadecene 


eloosene 


(keal mol) 


” 


(keal, mol) (keal mol) 


2-74 244 2-15 
2-50 2-38 
2-80 2-51 
2-06 2-72 
S11 2-87 
313 
319 


true of the alkylbenzenes, alkyleyclohexanes, 


alkyleyclopentanes and siloxanes [3, In 
general it may perhaps be said that E, and E, 
decrease with increasing depolymerization, irres- 
pective of whether this is brought about by 


chemical or physical means (temperature), 


APPLICATION OF CORNELISSEN AND WATERMAN’S 


VISCOSITY TEMPERATURE ForRMULA TO 


GLASS 


It has been investigated whether this formula 
can be applied to glass [7]. 


Table 3. 


Generally the formula is used in the form 
log 4 log » p A/T* +B 


The density of glass varies with temperature to 
a much lesser extent than its viscosity does, and 
therefore the factor 1 p has but little influence 
on the variation in viscosity with temperature. 
Hence log » can be used instead of logy. This 
will be clear from Table 3, which contains one or 
two examples extracted from the article by 
Suarrsis et al. [8]. 


Composition of glasses (in percentages by weight) 


Li, O SiO, 


20-0 


t ¢ low log 


2-200 


205 


0-004 


O-BS7 
0-338 


4 
1 1-82 1-58 
1 210 1-4 1-53 
| 2-69 2-50 2-25 1-08 
1 
3-65 3-43 
376 3-58 
Na, O K, O | ee ee | log p 
1101 
| | 1 | 
1-201 
23-9 76-1 900 3-523 3-166 
1400 1-396 1-058 
: 2-127 2-108 
17-8 $2-2 1000 2-933 2-504 
1301 1-430 1-104 
1-503 1490 0-013 
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Table 4. Glass of composition 3 SiO, - Na,O* log » 


Srevets and H. WarrerMan 


3-528 wea Te O-4747 


values 
Temperatur 
( K) Observed 
1574-6 
1551-2 
515-2 
1472 
1455 


1417 


‘ 


| 


Calculated 


1455 


> 5 


> 


i 


The values were obtammed 


Viscosity 


formula log » 4 B has 


composition 


In Table 4 the 
been apphed to 
Na,O nm a 
10? and 10* p. 

The 


silicate glass in a range of viscositi s between 107 and 


a glass having the 


ofl 


range viscosities between 


results of applying the sar formula to a 


10'* p are shown in Table 5. It should be observed 
that these Tables are only provided by way ol 
example, and that many more of the same kind 
could be given. 

The question now arises as to why formula (4) 
has a better degree of validity tha: formulae (1) 


and (3) have. When we were dealing ith hvdro- 


from the | 


boratories of Corning Glass Works 


of E 


decreases with increasing depolymerization and 


carbons we concluded that the value 


re spective ly mmereases with mecreasing poly ine riza 


Wi 


tion. shall now comsice r two Cases, 


(a) hermal depolymerization 
/ 


as its temp rature 


falls off 


ad to 


that the viscosity ol class 


fact 
increases must be 
the breakdown of the network (depolymerization 
It may be expected, on the analogy of this, that the 


of ky 


kor olass 


value will decrease at higher t rap ratures, 


then, the quantity KE ois a function 


of temperature, decreasing as the temperature 


: 

6 1918 

) 

aes 
0.2 

2 2 

o-4 

Os 

oo 

ot 

+ vol 

ad 

1200 26 

2s 
1523 
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Table 5. ¢ omposition of glass in percentages hy weight : 
SiO,: 672 Na,O OS KO 9-7 
CaO MyO: 2 A1,0,: 10 
BaQ: 2 ZnO 


38-10% 


Discre pancy 


Observed* Percentage 


SUS 
OLS 
114s 
154s 
Os 


libs 


*Lhe viscosity values were obtained from the Development Centre, Glass Division. N. \ 


Philips’ Gloeilampenfabrieken, 
The accuracy of the log Values ts Oog 


mereases. If its ck pendence on temperature ts as it does in the hydrocarbons). In other words 


expre ssed by the re lationship k 250 R. a decrease nh k accompanies an increase in the 
d (logy) d(1 T) ( T’', where C and a a number of ions modifving the network. 
constants, the formula of CORNELISSEN and \ similarity 


between chemical and thermal 
WaTrerMAN will be arrived at 


depolymerization is evident in the infra-red 


spectrum of glass: when the temperature is 
(h) Che mical de polymerization 


raised or the proportion of ions modifying the 
Not only does E depend on temperature, as network | 


s caused to increase, an analogous 


just noted, but we are faced with a phenomenon change in the intensity curves of the infra-red 


whereby it decreases with the ce vree of coherence absorption is observed [9, 101. 
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Heat transfer between a fluidized bed and a vertical tube 
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Abstract \l ‘ n the literature on heat transfer between a bed tMuidized bry 
im of a Tit vial sitioned tube are represented by two different correlations 
fine and light ! les where the viscous forces on the particles are predominant 

vtween the thu i bed and a flowing gas or liquid is assumed, and another for 

user and heavier particles where inertia effects prevail. The distinction between 

ilready bee le in a previous article on heat transfer t horizontal 

ituthonms ‘ i bilitw of the correlations are discussed and measures to 

<imizing th t transfer under given conditions are considered. The correlations 

betcanne vat te of the average coeflicient of heat transfer bet ween a thiidized 


madd ane vertical tubes which occupies a pre . of the cross-section of the bed 


Resume wuteur démontre que toutes les données qui se trouvent dans la litterature con 
chaleur ent in lit fuidisé au moven d'un gaz et un tube place dans lane 

representées par deux correlations quis ippliquent a deux inters illes differents 

illes, Dun « ti des grains fins et levers ou precdonmiune leffet des forces visq 

ueuses sur les grains et ot le lit f lise est considere comune nalogue a un courant de gaz ou cle 
lo autre intervalle des grains plus gros et plus lourds ou prey ilent les effets 
un autre ur Véchange de chaleur avee un tube horizontal, ill a cleja 
été publié en quoi ces deux int lies sont différents. Loauteur discute les limites d appieation 
des correlation conmsidere les esures nécessaires pour obtenir un éechange de chaleur mani 
run Sous conditions donnees Les corrélations aident aussi a évaluer le coeflicient moven 


dechange « chaleur entre un lit fluidisé et une batterie de tubes verticaux qui occupe une 


larue partie de la section trans 


Zusammenfassung \ile in der Literatur verfugbaren Werte fir den Warmetibergang zwischen 


einem Fliessbett. das durch en Giasstrom unterhalten wird und eimem senkrecht stehenden 


Rohr. werden durch zwei verschiece tezichunygen dargestellt die eine gilt fiir das Gebiet det 
feinen und leichten Kornungen die Adhdsionskriifte bei den Teilehen vorherrschen und ei 
Analogie zwischen dem Fliessbett und strOmenden Gasen odet Flissivkeiten angenommen wird 
withrend die andere das Gelbiet iberen und schwereren Kornungen, wo Tragheitswirkungen 
vorherrschen, erfasst. Der Unter tied zwischen diesen beiden Bereichen wurde schon tn ecimert 
friiheren Arbeit tiber den Wir vrgang in einem horizontalen Rohr darge stellt Dic Grenzen 
der Anwendbarkeit dieser Bezie! ven werden diskutiert und Massnahmen zur Erhohung des 
Wiirmetibergangs unter bestimmten Bedingungen erértert. Die Bezichungen konnen ebenfalls zur 
Berechnung eines Durchschnittswertes fiir den Warmetibergang zwischen emem Fliessbett und 
einem Biindel vertikaler Rohre diet das einen grossen Teil des Bet tquerse hnoitts 

1. Previous Work wrru Horizons Tun The one range is that of fine and light particles. 


IN eARLIER publications, a comprehensive series It is characterized by 


of measurements of heat transfer from Huidized GD 
beds to horizontal single tubes has been re ported . 
1, 2). A simple correlation of the results of the re 

measurements made before 1951 was first given \ new correlation for this range was derived 
»). but it had to be revised when more data’ from an assumed analogy between the thuidized 
became available [1] and it proved useful to bed and a tlowing gas or liquid [1 It contains 


distinguish between two ranges. the following dimensionless yroups 
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Heat hetween a 


transfer 


AD, 


characte rized by 


of coarse and heavy partn les Is 


The correlation for this l contains 
the following dimensionless CTOUPS 
AD, GD p (ji 
and 


In the intermediat 


range 


GD 
2550 


pie 


2O50 


it was recommended that the average of the heat 
transfer ciM tlic ments predic ted thr two corre la 


tions be used. 


Re mark 
The 


place not only when D, or 


the take 


p, IS varied, but also 


transition between ranges can 


when gas properties like G, p, or are varied. 


Theretore is it incorrect to Spe ak of ranges of fine 
and light materials on the one hand and of coars« 
materials on the other. But on 


and heavy may 


uphold this convenient fiction as long as one 
the effect of the 


and other vas properties. 


bears in minds mass velocity 


2. Previous Work wiru Vertica. Tures 


Table 1 all the 


heat transfer between beds Hhuidized by means of 


enumerates investigations into 


gases and vertically inserted evlinders. There is 
no uniformity in the methods used by the authors 
to determine the average particle diameter. It is 
impossible to account for these differences sincé 
the particle size distribution is not given. The 
appendix gives additional details 

A simple correlation of the ¢ xperimental results 
up to 1951 [3, 4, 6, 8 


where [8]. 


has been prese nted else- 


When, however, more and more experimental 


data became available 15, 7, 9, 10, 11, 12 
of these [7, 9, 11] 


old correlation [8]. 


SOTTI 


necessitated a revision of the 


fluidized 


Erperiments with avial cylinders 


Table 


bed 


Va 


diameter 


Diameter 
culinde 


diamete 


ticle 


Reference 


ina 


verth 


tule 


Variot 


ZI O NTE 


ir 


Ou 
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i 


beads 


co) 


ALO) 


16-7 
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Sand 


100 


7 0188 
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CorRELATION FormMuLar ror Heat TRANSFER 
rrom A Fiucipizep Bep ro Verticat 
It proved useful to distinguish between the 
same two ranges as in the case of horizontal tubes. 
\ new correlation for the range of fine and light 
materials can be derived from the assume d analogy 
with a gas or liquid flowing through the annular 
space between the evlindrical outer wall and the 
axial tube. 
Watcer [13] has published a critical survey of 
} 


the literature on heat transfer between the walls 
of an annular space and a gas or liquid flowing 
through the Hk 


that coetlicrent 


annulus, 


tha 


turbulentl,y arrives 
at the of heat 


transfer at the inner wall can be represented by a 


com lusion 


simple relation between the Nusselt group, the 


Reynolds group, the Prandtl group of the fluid, 


diameters of the outer and 


inner evlinders, and the 
ties of the fluid at the bulk temperature and at 


the ratio between the 


ratio between the 


Visco 


the temperature of the inner wall respectively 


\. DENBERG 


the difference between the diameters of the two 

evlinders has to be taken the 

length in the Nusselt and Reynolds groups. 
Conseq ue ntly the assumed analogy leads to an 


to the 


characteristi 


as 


attempt 
heat 


and light particles and Ve rtically imserted « \ lirncte rs 


represent exp rimental results 


on transfer between fluidized beds of fine 


by n re lat won by tween 


h (D, 


dD, G(D, D,) 


The 


these 


the first three of 
the 
ponding correlation for horizontal tubes (Section 
l the of the factor (1 the 
Revnolds group. It be omitted 
the porosity is not given in most of the publications 
listed in Table 1 


motion and consequent 


there 


mam difference between 


and those occurring corres 


Is absence € ih 


has to because 
In view oft thre violent particl 
of 
the 


at an other than the bed te rature 


uniformity tempera 


im no reason to introduce viscosity 


ture 


of the 


Correlation of results of experiments in Tables 1 and 2 (see also Fig. 2) 
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The investigations enumerated in Table curves thus obtained ; they include the transition 
comprise 110 series of measurements. In twenty rang 
two of these the experimental conditions wer GD, 
2050 
entirely within the range 
where the correlation for fine and part les 
"050 
has to be used in combination with the one for 


coarse and heavy materials. The broken lin 
which we qualified as the range of tine and light : 


D, — D,)p,\" 
conditions were at least partly within the above i] (5) Cu on | 


articles, In another twenty-four series the 
}, ), 


range. It was found that a fair representation 
of the forty-six series or parts of series within the where 
range can be given by plotting thy composite Sand B 
group 
h(D, — D,) (D,\"" (k\'" G(D, — D, 
nie 


asa tion of thre he vnolds group 2 ind RB whe n 


7 G(D D,) 
G(D, D,) 0-287 108 


Figs. 1 and 2 show the forty-six experimental has been drawn through the set of exp rimental 


results of experiments in Tables 1 and 
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curves: the co-ordinates of any pont of these 
curves differ by no more than a factor 2-2 from 
those of the nearest point of the broken line. 


\ correlation for the range ot coarse and heavy 


particles can be obtained by adding two dimen- 


sionless yroups to the four oc urring in the corre 
lation of Figs. 1 and 2. namely 
am 
g dD, dD, 
It was found that a fair representation of the 
eightyv-eight series or parts of series of measure 
ments within this range and the transition range 


can be given by plotting 


h(D, — D)(D, D, 


as a function of 


G(D, 


The latter composite group may be written in 


the simpler form 


and is thus found to be equal to the product of 
two other well-known groups, viz. the reciprocal 
of the group D, (D, D),) introduced above, 
and the square root of the 
D,, p* gil. 14 

Figs. 3 and 4 show the eightv-eight exper 


Froude group 


mental curves. The broken line 


h(D, D,) dD, D, 


where 


?-O when 


O105 and F 
G(D, 


p 


Correlation of results of experiments in Tables 1 and 2 (see also Fig. 4) 
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hig Correlation of results 


and E 240 and FF O-10 when 


G(D, 


Avain. 


curves 


has been drawn through the set of curves. 


the co-ordinates of any point of the “i 


differ by no more than a factor 2-2 from those of 


the nearest point of the broken line. 
The 


curves 


experimental 
Figs. 1-4 
tubes [1 


deviations between the 
and the 
are larger than 
But then the latter 


experiments only, 


correlation lines in 


those for horizontal 
were derived from our own 
there are no other data 


available in the literature. whereas the former are 


based on experiments by a great 


authors as specified in Table 1. There are wide 
differences between their CX rimental set Ups ; 
the main dimensions are given in Table 1. It is 
that the lll, 15 


10] of the axial heating or cooling 


known dimensions and con 


figuration [4, 


number of 


experunents in Tables 


method of introducing 


affect the 


evlinders, as well as thr 


the fluidizing 12). vreatly experi- 


mental results, The greater part of the differences 
Figs. ] J 


may be ascribed to these influences. which cannot 


between the rimental curves in 


be described quantitatively in detail as is th 
ease for the variables contamed in the correla 
tions, because the number of systematic experi- 
ments concerned ts too small. 

In their atte mpt to establish a general correla- 
WENDER 
restricted number of the experimental results of 
Table 1 [3. 4. 7. 


between expert- 


tion and Cooper 16 included a 


some of the authors listed in 
11 The 
mental data and their correlation line are 
in Figs. 1-4. 

When the mass ve locity G is increased at con- 
bed 


rapid 


deviations these 


much 


larger than those 


stant gas density and temperature, the 


correlations predict a increase in fA in 


Part I of Figs. 1 and 2 and in Part III of Figs. 3 
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Il. A. 


and 4. The values of G involved ther will not 
be much greater than the minimum value at which 
In Part Il of Figs. 1 and 


” there will be a moderate increase in / Owing 


fluidization can occur 


to the higher values of G which corre spond to this 
region, fluidization may be more pronounced 
than in Part I, but the bed is still assumed to 
gas or liquid. In Part lV 
a shght 
takes 


paths 


behave like a flowing 
ol Figs 3 and however, there will 
fluidization 


decrease in hk. Here intense 


place with high bed porosities 


from those 


of the partie les aifter essentially 


ras number of experime! curves 


obtained with mecreasing clon clearly 


show the transition from a rapid inereas h 
nerease to a shi 
thers do not cover more tl 
rtwo [3.11] of these stag to the 
ittains un Valu 


when G is mereased t it the 


transition to Part IV of by 
place At lower values of G the cond 
represented either by Part Ill or by some part 


correlations alw i\ 


takes 


ol Figs and 


e mass veloc t which 


maximum Value follows 


the same time ondition 


D D 


has to satistied Since 
ecorretations ss imect to an uncer 
be advised | ec a mass velocity 
value corresponding to the nun 


h in order to make certam | transition 


to Part IV will alread nave | Ther 


the result will by » coelhoent transter 


0-03 times the maximum value 


equal to 
according too the correlation : thy erence is 


unimportant in w of the uncert 
application of the correlations. 

When the diameter D, of the ax! evlinder ts 
velocity, bed mperature, 


bed diane r oar kept 


increased while mas 
particle propertic 


constant, the correlations pre cict 


h in Parts Land ILL, the same in Part IL as long 
as D, dD, »-68. and also in Part IV as long as 
D, < Dy, 3-30. 

Since Table 1 shows that in all « xperiments on 
which the correlations are based the ratio 
DD, [3 


them when the ratio is as low as 2-68 or 3°30 


it is not pe rmitted to apply 


anyhow, 
for values of D, D, such that th 


thev always pre clict 


Therefore, 
correlations are applicable, 
decrease in hk. It follows that the use of tubes of 
small diameter is always to be recommended 
when it is desired to obtain high heat-transfer 
coethcrents. 


Similar statements about the effects of other 


yy rtinent variables on 4 can easily be derived 


from the correlation formulae. In practical 
applications it will seldom be possible to vary a 
quantity such as G, D, and D, whale keeping 
the others constant If certain relations between 
these quantities are given it ts, of Course possi bole 
to derive from the correlations the combined 
influence on /. An example will be worked out 


my the 


Rance or Vauiprry or THE CORRELATIONS 


In view of the derivation of the correlation for 
fine and light particles, it 1s permissible to apply 
it when the relevant dimensionless yroups ar 
within or nearly within the range cove red by the 
experiments, Therefore it ts recommended that 


it’ bee only in such cases where 


the Reynolds group GD. smaller than 
2050 | 1 
group 


and 5 


the Revnolds 


between O-1 
the Prandtl vroup ¢ ko iw between O06 and 
13: 


the ratw is between 4-0 and 


In the range ol coarse and heavy particles the 


have been introduced only as a means of arriving 


at a useful correlation of the exp rimental data. 
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It follows that its application can be recommended 


only when ; 


the Reynolds group GD, p,/ pj > 2550 ; 


all the variables are within or nearly within the 
range covered by the experiments as specitied 


in the last column of Table 3 in the appendix. 


For values of the Reynolds group GD, Py py 
between 2050 and 2550 it is recommended that 
the average of the values given by the two 
correlations be considered as a good prediction 
of the heat-transfer coetlicient h. 

The correlations in Figs. 1-4 apply to heat 
transfer between fluidized beds and axially 
positioned cylinders only. Coeflicients of heat 
transfer to vertical tubes in eccentric locations 
are about 1-6--1-8 times as high, when the distance 
between tube and axis is 0-35-0-71 times the 
bed radius [8]. Comparative measurements of 
coetlicients of heat transfer at the wall of the bed 
and in the axis [4, 10] have shown that, depending 
on conditions of fluidization, the former may be 
either higher or lower than the latter. but both 
are lower than those in intermediate locations. 

It is Impossible to make a good estimate of 
the magnitude of these effects because the 
differences between the results of various experi- 
ments on heat transfer between fluidized beds 
and the wall [16, 17] are even larger than those 
for axial cylinders, In any case the influence of 
the eccentricity of an inserted tube seems to be 
smaller than the deviations between some experi- 
mental curves and the correlation lines in Figs. 14. 
Therefore, when a reasonable estimate of the 
average coetlicient of heat transfer to a bank of 
vertical tubes which vecupies a great part of the 
cross-section of the bed is to be made. it may be 
recommended that the correlations as given in 
Figs. 1-4 be applied. 

Finally, it might be interesting to check how 
far the correlation established for cases where 
gas is the fluidizing agent is also valid if this 
function is performed by a liquid. Heat-transfer 
measurements in this type of system have been 
reported by Wessex and Marpus [18], who used 
water as the fluidizing medium. However. in their 
case the Reynolds group G(D, — D,) p,/pyu was 


smaller than 0-02 10° so that their data are 
outside the scope of our considerations. The 
production to the left of the straight line in 
Part IT of Figs. 1 and 2 happens to be a reasonable 


representation of their experimental results. 


5. APPENDIX 
§.1 Experimental data and conditions 

Table 2 specifies the experimental data pertain- 
ing to each curve in Figs. 1-4. Those experiments 
of Mickiey and in which the bed 
porosity was higher than 0-70 have not been 
included, since they are outside the scope of 
the present investigation [8]. This limitation 
does not affect any of our conclusions. OLIN and 
Drax [7] paid special attention to the influence 
of fouling of the particles on heat transfer. 
Only their results obtained with clean particles 
have been considered. The data obtained by 
Wicke and Frrrine [9] with activated charcoal 
particles have not been included since the particle 
density is not reported. They state that in three 
other cases fluidization was disturbed by slugging 
phenomena. These cases have also been excluded, 
although the results obtained with SiC particles 
iluidized by hydrogen do agree with the correla- 
tion. Of course, data obtained at mass velocities 
so small that fluidization could not occur [3, 5, 9, 
12] have not been considered either. 

Electrical heating elements were inserted in the 
axis of the bed by the greater part of the authors ; 
in two investigations [6, 8] a water-cooled vertical 
tube was used instead. 

Table 3 gives the extreme values of the most 
essential variables as applied in the experiments. 
Data for fine and light particles are given 
separately from those for coarse and heavy 
particles. 


5.2 Example of application of the formulae 

Suppose that, for a certain annular bed of sand 
which is to be fluidized the air mass flow rate be 
®,,, and that the pressure, the temperature and 
the degree of fluidization be given. 

Then the mean diameter D, of the sand particles 
and the diameters D, and D, of the cylindrical 
walls required for maximum coefficient of heat 
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transfer between the bed and th 


be deduced using our correlations. 


inner wall can 


It follows from the data given for various kinds 
of sand [1, of the 


product of mass velocity and kinematic viscosity 


2) that the minimum valu 


at which fluidization of the sand could be observed. 
IS approximately proportional to D.*?, and that 
the proportionality 


0-26 kg m*? sec?. 


constant > IS 


equal to 


Therefore the specified degree 


of thuidization is conditioned by : 


where 


« iS a given dimensionless number which is 


greater than unity. The given mass rate of flow is 


expressed by 7 


D?) G 


m 


(D2 
2 


Elimination of the mass velocity G gives: 


ps 
TT k 


p 


When the expressions found for G and D, are 


substituted in the correlation formulae, the result 


for the range of fine and light particles reads : 


h(D, — D,) | ) 1! p, 
D, la (D, D,) py} 
or 
H D,\'* ! 
D, p,\n,) (D, D,)" 


For the range of coarse and heavy particles one 


has 
h(D, D,)*3(D, k ) 

h dD, cy 

(D,* D,*) p} 
or 

A. 1/3 
h=J(D, ve M2) (D, + D,)*” 


D, 


Here H and J are constant as long as D, and 


D, are the only variables considered. 


If. in the first place, the effect of an increase in 


D, at a constant value of D, is considered, the 
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formula for fine and light particles is found always 
to predict a decrease in h. 

For coarse and heavy particles an increase in 
h is always obtained in Part III of the correlation 
graph where F 2-0, and a decrease is predicted 
in Part IV where F 0-10 as long as 
D, < 69-4 D,; this mathematical condition is 
always fulfilled when it is permitted to apply the 
correlation, since the expermmental basis covers 
the range 4-0 D, dD, 30 D, only. 

The increase in dD, is, of course, coupled with a 
decrease G which is proportional to 
1 (De D,*), and a decrease in D, which is 
proportional to (D,? 


the conditions remain within the range of fine and 


As long as 


light particles, it is desirable to decrease the value 
of D, as far as possible, i.e. until the transition 


range is reached, At this limit we have 


5/3 
GD, Ps | Ps 2050 
py (D,? D,*) p} 
If there 
«B(D, D,) 
L070 
g} 2 


so that Part IIL of Figs. 3 and 4 comes into play, 
a further decrease in D, would have a negative 
effect, the maximum value of A having already 


been attained. In Part IV where 


dD, 


gh? 


) 
D,) - 1070 


higher values of h can still be obtained by means 


of a further decrease in D,, until the maximum 


ordinate in Figs. 8 and 4 is reached: Le. 
«B (D, D,) = 
1070 


Summarizing, the most favourable value of D, 
at constant D, is determined by either 


®,, 


Ps 2050 
(D2 D,?) p} 
— Di) — 
pg? 


a 

. 
te 


H. A. Va 


5 . 


«B(D, 


git 


D,) 


1070 


In practice, one will solve each of the equations 


| ?,, Pe 
(D,* D,?) 


1070 


for D, ; the smaller of the two values thus found 


is the correct solution. 

In the second place, the effect of an increase In 
of D, can bye 
and light part 
D, 1-61 im Part I 


and as 


D, ata constant value considered. 


The formula for fine cles predicts 
a decrease in h as long as D, 
of the correlation graph where B 3-4, 
D, 3°29 in Part Il where B O-44: 


these mathematical conditions are fulfilled when 


long as dD, 


it is pe rmitted to apply the correlation, Le. when 
+0 D, D, < 25 D,. 

For coarse and heavy particles a decrease h 
Is always obtained in Part ILL of the correlation 
graph where F 2. and also in Part [IV wher 
F 0-10 as long as D, D, 2-96 ; 
condition is again fulfilled when it is pr rmitted 
30 D,. 


this 


to apply the correlation, 4-0 D, dD, 
The increase in D, involves an increase in G, and 
an imecrease In D,,. In all cases it is desirable to 
decrease D, as far as possible ; a lower limit may 
be set for mechanical reasons or by the require 
ment that the inside of the vertical tube should be 
accessible for cleaning or maintenance. 

The answer to the question of how to obtain the 
maximum value of h is, therefore, to choose the 
diameter D, of the 
possible ; then the diameter D, of the outer wall 
can be found as outlined above. The corre spond- 
ing values of G and D,, can easily be dete rmined, 
The solution thus obtained is valid only when 
< 80, and when D,, ts within or 
nearly within the region covered by the experi- 
} has been derived, 


inner evlinder as small as 


ments from which the value of 
in other words when the mean parti le diameter 


is between 80 and 500 jx. 


Nume al rample 


Let the following numerical data apply : 


30 
0-0408 
0-67 

1030 
O-OOL0 


OO, 


while the minimum allowable tube diameter dD, 
is 0-020 m. The 
and 4 is then reached for D, 
these values of and 


ordinate in Figs. 3 


00-0844 m; for 


maximum 


®,, 
(D? D,?) 


1070 


The corre sponding values for mass velocity and 
mean particl diameter are G 0-189 kg m?* sec 
and D, B04 10 Since D, Dy 
is between 4-0 and 30, and the particle diameter 


is between SO and 500 ja. this solution is correct. 
In this case # follows directly from the maximum 


value of the ordinate in Figs. 3 and 4: 


hA(D, dD, 


D, 


whence A 668 Jo see 

In the preceding discussion of the formulae 
it has been recommended that the value of the 
abscissa in Figs. 3 and 4 be chosen twice as high 
as the value corresponding with the maximum 
ordinate ; the reasons were the uncertainty tn the 
actual position of the maximum and the fact that 
the ensuing decrease in the caleulated value of h 
was as low as 7 per cent. Although in the present 
example the decrease im h would be 8-5 per cent 
because of variations in D, and D, occurring 
in the ordinate of the graph, this might induce one 
to choose here also «8 (D, D,) pg? 2140 
instead of 1070. In that case, however, D, would 
become so large that the conditions would already 


be in the range of fine and light particles, It is 
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Heat transfer between a fluidized 


therefore certainly not advisable to increase dD, 
beyond the value at which the range of fine and 


bed and a vertical tube 


NOTATION 


1 proportionality constant 
light particles is entered, this value being deter- B = exponent 
mined by ( specific heat of particles at bed temperature 
J/kg 
| Ps 050 D, mean volume-surface particle diameter [19] m 
(D2 D,2) p| dD, outer diameter of diametrical tube m 
D, outer diameter of axial evlinder m 
Thus dD, 0-124 m, «8B ( dD, D,) 730, D, diameter of fluidized bed m 
G 0-086 ke /m? see, D, 182 10-8 m. For proportionality constant 
these values of D, and D, the heat-transfer ” = exponent 
mass velocity of fluidizing gas kg m® se« 
coellicient is 623 J sec according to the 
‘ aes g acceleration of gravity m sec* 
correlation for coarse and heavy particles, and I — proportionality constant JmB-1! ‘sec 
647 J, m*? sec °C according to that for fine and h coetlicient of heat transfer between bed 
light particles. The recommended estimate in the and tube J/m? see “¢ 
yroportio consta F+16/9 
transition range is the average of these two J = proportionality constant J/m wes 
d thermal conductivity of gas at bed 
values; it amounts to 635 J m* see °C and is 
temperature J see m 
lower than the maximum value calculated above proportionality constant kg /m3!2 see? 
by 5 per cent only. « — fraction of voids 
It may be concluded that, in view of the x — measure of the intensity of fluidization 
uncertainty in the correlations, it is indeed # — dynamic viscosity of gas at bed temperature 
advisable to choose D, 0-124m,. G 0-086 ko nae 
. = density of gas at bed temperature kg m°? 
m® sec, D, = 182 10°* m instead of the values p, density of particles 
calculated above. ®., mass rate of flow of fluidizing gas ky sec 
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The motion of a rigid sphere in a frictionless cylinder 
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Abstract This investigation was undertaken in order to devel pa opohne re and evlinder model 
which would serve as a basis for further theoretical studies of assemblages of parti les Stokes 


equations of viscous flow were assumed to apply to the steady motion of a rigid sphere slong the 


axis of an infinitely long evlinder. Solution of the problem is carried out under the postulate 


that the fluid shear on the cvlinder walls is everywhere equal to zer 


An exact solution is obtained in the form of an infinite set of linear simultaneous equations 


for the coeflicients in the Stokes stream function 
The above problem was solved for a sphere to evlinder diameter ratio of O-1 to ¢ 7. Results 
are presented in the form ofa plot of the stre un lines from the solution of the Stokes’ approximation 


in the viscous flow region Applications t flow relative to assemblages of particles ts dise ussed 


Résumé Cette ctuct treprise dans le but de «te velopper un motele ce sphe res et de 
mibl os ch particules Les 


eviindres, pouvant servir 1 une étude théorique des ass 
equations de Stokes relatives 4 des ccoulements visqueux, sont supposees sappliquer au motive 
ment permanent d'une sphere le long axe dun eviindre infiniment long. La solution 
du probléme est effectuée en supposant que le cisaillement d'un thuide sur les parois du « vlincre 
est partout égal a zero 


solution exacte est olbbtenuc sous forme dune serie equations lneaires simultanees 


pour les coefficients de la fonction d écoulement de Stokes 

Le probléme ci-dessus a été resolu pour un rapport du diametre de la sphere au diametre 
du evlindre de 0.1 4 0.7. Les résultats sont preséntés sous forme dun diagramme des lignes cle 
écoulement provenant de la solutrar ippreximation de Stokes dans ha region che Pecoulement 


VisqueuN., Les applications a un écoulement coneernant les assemblages de particules sont 


discuteées 


Zusammenfassung ies niersuchung wurele Entwicklung eines Kugel- und 
Zviindermodelis unternommen ins «als Grundlage fiir weitere theoretische Studien tibet 
Teilchenansammilungen dienen soll. Die Stokes’ schen Gleichungen des viskosen Fliessens wurden 


ingenommmen fur di lei hférmige Bewegung einer starren Kugel entlang det 


als anwendbar 
Achse eines unendlich langen Z5 lers. Die Losung des Problems wird unter der \V oraussetzung 
durchgefiihrt, dass die Fliess-Beanspruchung det Zviinderwiinde tiberall gleich Null ist 

Kine exakte Lésung wird erhalten in Form einer unendlichen Reithe linearer simultanet 
Gleichungen fiir die WKoeflizient les Stokes schen Stromungsyeset zes 

Das erwihnte Problem wurde fiir cin Durchmesserverhiltnis von Kugel und Zvilinder von 
0.1 bis 0.7 erhalten. Die Ergebnisse sind in Form cines Diagramms der Stromlinien der Stokes 
sehen Niherungslésung in zihen Bereich angegeben. Anwendungen aul die So romung re lativ zu 


Teilchenansammilungen werden diskutiert 


1. INrrRopUCcTION several investigations in the past. 


investigated the effect of a rigid sphe re moving 


Tue DEVELOPMENT of a “ model cell arrange 


ment which could be used as a basis for theoretical through a viscous fluid along the axis of a tube 


studies involving pressure drop, tuidizing velocity, employing the Stokes -Oseen resistance law 


ete. m packed beds has been the subject) of Wakiya [2] treated the case of stationary 


: 
ig 
vol 
1050 
* 
i 
: 


sphere in the flow of a viscous fluid in Poiseuille 
flow. Both investigators gave an approximate 
expression of the and 


sphere drag. 
[3] arrived at a solution for the case of a 
moving sphere in the path of a moving viscous 
fluid by the use of a method of * reflections.” 
In this case the velocity field was decomposed 
into its components, each of which satisfied the 
wall boundary conditions exactly and the sphe re 
boundary conditions approximately, 

A recent paper by Richarpson and Zaki [4] 
proposed a cell model derived from the assumption 
that the sphere was surrounded by a hexagonal 
envelope of fluid. They were then able to idealize 
this assumption for the case of a cylindrical 
envelope and thus give an expression for the 
approximate drag of the sphere. 
method of 


HawermMan [5, 6] 


presented a 


solution for the case of a stationary or moving 
sphe re located axially na cy line r through which 
a viscous fluid flowed, Cylindrical and spherical 
Stokes’ 


flow were used to satisfy the boundary conditions 


co-ordinate solutions to equations of 
on the evlinder and sphe re respectively. The 
cylindrical boundary in this case was a physical 
wall and hence was not frictionless. Kawacut! 
[7] investigated a similar case, using the method 
of reflections” similar to Harrre: and Byrxy 
[8]. An approximate expression for the drag of a 
sphere in a frictionless evlinder was obtained, 
applicable up to a sphere to cylinder diameter 
ratio of 0-2-0°3. 

The problem 


discussed in this paper was 


solved using a modification of HaperMman’s 
approach. The results are then exactly applicable 
throughout the entire range of sphere to evlinder 
ratios, unlike solutions derived in the past where 


methods other than HanerMan’s were employed. 


Morion OF A SPHERE IN A Fricrion.ess CYLINDER 


The thiid system investigated here is assumed 
to have such characteristics that Stokes’ equations 
of motion apply. The Stokes’ stream function 
expressed in terms of cylindrical co-ordinates is 
used to satisfy the boundary conditions on the 
evlinder walls and at infinity. This expression is 
then expanded, using both co-ordinate systems 


for convenience and compared with the stream 
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function 


expressed directly in spherical co- 


ordinates. To enable the above comparison to 
be made, part of the latter expression is trans- 
formed into the cylindrical co-ordinate system 
where necessary. Comparison of the terms then 
yields a relationship between the constants. The 
boundary conditions on the sphere vield rela- 
tionships between the constants in the spherical 
co-ordinate solution. Substituting the previous 
relationships into the relationships obtained from 
the sphere boundary conditions then yields an 
infinite set of linear simultaneous equations for 


the solution of the constants. 


Theoretical deve lopme nt 


The co-ordinate origin is taken at the centre of 
the sphere and the cylinder is assumed to be 
moving at a constant ve locity U in the negative 


x-direction (Fig. 1). 


be ve ty rutorm fluid 


Fig 1. Co-ordinate svstem. 


The boundary conditions are : 


On the surface of the sphere (r R): v, 0, 
Vg = 0. 
At infinity (4 0): v, U, vp = 0. 
At the evlinder walls (p = b): = = 0, op = @. 
) 
dp 


The Stokes’ stream function expressed in terms of 


eviindrical co-ordinates ts [5]: 
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s Substituting for the A, and A, integrals and 
db | [p K, (ap) f, (a) + p® Kg (ap) Fy (a) expanding the remainder of equation (6) In a 
: Taylor series we obtain: 


p I, (ap) (a) 


77 


p (ap) G, (a) ] cos axvda (1) rsin® + ’2 sin* @ (— 8cos?@ +1) 
> 


2 
velocity components are : ' 
1 1 sin? @ { 15 cos? + 9)cos@ 


pop p oa r? sin 
1. 
sin? #(945cos' @— 630 cos? + 45) 
rsn@odr 
Satisfaction of the boundary conditions on the 
evlinder walls vields the following two equations 
aK, (ab) f(a) [2K , (ab) abK ,(ab) | F(a) 
al ,(ab)g,(a) [2/, (ab) abl |G,(a) 0 
where : 


K, (ab) f,(a) 4 bK,(ab) F,, (a) I, (ab) g, (a) 
bl, (ab) G, (a (2) [g, (a) a 2 G, a" da 


We can therefore solve for g,(a) and G, (a) in 0 


x 


G, a” da 


terms of f, (a) and F, (a): 


kK, (ab) 


(a) 0 
I, (ab) 


(@) 

The stream function obtained directly from 

Bi a (7, (ab) P 

(r, @) c 


Letting 


K, (ab) 


we obtain: Expansion of the above equation and term by 


g, (a) S, F, (a) S,f, (a) term comparison with equation (7) vields ; 
G, (a) S, Fy (a) (5) 


Expansion of the cylindrical co-ordinate stream 
function yields : 


. 


(a2, p) ph, (ap) [ay + a,a +a,a° 


0 COs aa da 


K, (ap) [bo t+ b,a + b,a* .| 2 (2n 


cos axrda 4 2 p? we obtain also: 


x 


[ (ap) a, (a p?I,(ap) G, (a) | 
cos azda + 


4 
(7) 

VOL 

> 
: 
oe 
(4) +B ** +D (9) 
r 
1 
D, = by 
4 
x 
Q 
a4 
(6) 1 
2 2 5 
‘ 
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From A,, A,, C,, C, we get the infinite set of 


n 


linear simultaneous equations for bp, iveuas a 


S, (ab)" d (ab) (13) 
d / 
Ss," | S, (ab)” (14) 
= 


and from the boundary conditions on the sphere 


we have: 


Pr 
n=2,4 


Hence 
2n +1 l 2 l l 
A, B. D, (16) 
and : 


= Ss,” ! \™ 2, (18) 


ete. 


“3 [— 3-225 
h b. 
[ 0-82 1 — 4-755 A” 
R 5 
ay 
1-6125 A 
[ 
<3 9-75 A*] 0 


b, 
[— 0-17916 a5] 4 2 [62 — 1-82 


“3 [427] +... (20) 
RS 
ete, 

For the present purpose it was ascertained 
that the infinite series could be evaluated with 
suflicient accuracy to determine the drag coeffi- 
cient by using eight simultaneous equations as 
shown in equation (20), 

In order to test the applicability of the method 
numerical computations were carried out for a 
range of A from 0-1 to 0-7 using an I.B.M. 650 
computer. 
For example : 


for A tub 0O-land R 1, the constants are : 
bh, 0-57340 U 
b 0-00021622 U 


a, 019116 U 
0-0000380881 


The drag on the sphere has been shown to be 


(5). 


d b, (21) 


A, = (— 
2)! 4)! (2n 3) 
C.=-h. — 352A) + — — 2-60 4 
B 
ws 12) a, |o7 
Letting 
VOL, 
59 7 
240 A") + [82 — 8-225 A5] 4 
Ri? Re 
i A x 
2,4 
r 
~ 
+(n (n — 3) D, sq] (15) 
r 
c x 
R, (n 2)! (2n +1) 
b 
| R™ A (19) 
m=0,2.. 


J. Harrrt 


\ wall correction factor can be defined on the 
basis of the drag on a sphere im an infinite medium 


(i.e. D Or RU) as follows: 


k bo ve 
Gn Rl 


Thus we have 


Table Correction factors to Stokes’ drag 


S pie 
cyl 


diameter aliw (eq 


Solids com 


177 
O-O08 1-420 


2450 


125 
188 


In the KAWaGcutis [7 theoretical 


treatment is available for comparison with our 


dilute range 


result. He gives, 


Ak (23) 


164931 A O-TOGRO A 


For A Ol. value for 
K 1-108, which is in excellent agreement with 
Table 1. 


of the reflection technique in the dilute range. 


i 


equation (233) 


the value from indicating the validity 

No theoretical treatment or data are available 
for directly checking the value of AK in the more 
It is difficult to see how a 


frictionless tube could be set up expe rimentally 


concentrated reoion. 
result for a sphere m an 
evlindrical duct 
These data (9) not re ported in detail 


so Haperman’'s [6 
actual was checked expert 
mentally. 
here showed agreement to within several per cent 
of the theoretical result in the cas of A Or, 
for the drag on a sphere suspended in a viscous 
fluid. 


It was also thought to be of interest to establish 


the nature of the (hid flow pattern for this model, 


and P. A. As‘ 


since the rate of decay of velocity in an axial 
direction should represent the most extreme case 
encountered in sedimenting systems (1.e. the case 
for no spheres above or below the one we choose 
This 


numerically for the case of A 


to consider). evaluation was carried out 


0-6, using the 


values for the constants by, b,, a, and ay. 


Fricthioniess wall, sphere moving in stil fluid 


Fic. 2. Streamline plot ~ 


Thus we obtain : 


3 
1466 0 
0-367 U 
0-0486 U 
dD, 10-750 U 


D, = + 2-087 U 


To check the accuracy of the constant evaluations 
we find the velocity forward of the sphere at 
l: 


200 


) 
4 
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a 
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cos @ | B. 


(5 cos’ @ 3 cos @) 


+ D, 
cos @ 


0-95 95 per cent of which is suffi- 


ciently close. 


Substitution of the constants 


vields the streamline plot shown below (A O-6 
Of special interest is the axial velocity com 
ponent v,, forward of the sphere. For this case 


substitution in equation (24) vields the following : 


diameters 


Sp) ‘rT 


on he re 


Thus the velocity is negligibh distance over 


1-3 diameters from the sphere. 


Application 
The 


calculation of drag on an assemblage of particles 


development above is applicable to the 
based on the moc lofa eviindrical fluid nvelope 


surrounding each particle, if an appropriate 
assumption can be made relating the model to 
the concentration of solids in the 


The Rb A 


diameter to cylinder diameter. 


assemblage. 


value of only fixes the sphere 
In addition it is 
necessary in effect to establish the cvlinder length 
of a unit cell to obtain a relationship between A 


and the solids concentration 4 of the form: 


(26) 


ZAKI 
and by Kawacuti give values of k ranging from 
0-94 to the arrangement 
We have taken 4 l 
here on the basis that the effect of the shape of 


Assumptions made by Ricuarpson and 
1-03 depending on 


assumed in the assemblage. 


the outside boundary is not important. This same 
relationship would apply in the case of a spherical 
outer cell boundary. Values of 

toh 1 are tabulated in Table 1. 
A in this Table then the 


ricnce by a sedimenting sphe re ina SUSpt nmsion 


» corresponding 

The values of 
vive resistance 
of solids concentration as compared with the 
drag on a sphere in an infinite medium. Thus. 
concentration — of 


so that the 


for example, at a_ solids 
@ = 0-125, K — 3-991 
sphere in such as ass mblag would be predic ted 
equal to 3-991 (Stokes’ Law drag). 


words the spheres in such an assemblage would 


drag on each 
Or in othe r 


sediment at a rate equal to 1 (3-991) as fast as if 
they were settling in an infinite medium. 

It is interesting to compare the predictions of 
the sphe re-evlinder model employed in this 
paper with a similar mathematical treatment [8], 
which has been shown to be in vood agreement 
with a substantial amount of experimental data 
pre viously published. This latter model is based 


on the assumption that two concentric sphe res 


issemmblawe 
The 


particles in 


can serve as the model for a random 
Inner 


the 


of sph res moving relative to a fluid. 
sphe re comprises on of the 
consists of a 
This 


the 


mblage 


Hhuid 


and the outer sphe re 


envelope with a “free surface.” 


value for 
to Stokes’ Law: 


treatment gives the following 


correction factor A 


kK’ 


Values of A 


reasonably well with values of A predicted on the 


obtained from equation (27) agre« 


basis ¢ 


to Ab 


by the sphere evlinds r model, up 
0-216, where A is approximately Il per 
cent below the corresponding value of A’, 

Our results do not agree very well with thos 
derived theoretically by Ricnarpson and Zaki 
[4]. who obtained substantially smaller values 
Stokes’ factor A. Their 


theoretical treatment assumes that the spheres in 


for the law correction 


: T 4 1, +. r? D > 
r 
ays. 
rs (24) 
V OL ne2,4 
ll D (28) 
| | 
2 os oes 
346 1-30 
. 
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an assemblage will be lined up directly one above NOTATION 


the other and this may account for the smaller 
A constant 

resistance to flow thus predicted. These authors 
constants in evilindrical co-« 


also obtained experimental data on fluidization sobution 


which is in agreement with their predicted 
constants mn rical co-c rdimates 
relationships. Their results are discussed and 
‘ sole ‘ 
compared with other data in reference 8]. Review radiu ylinder 
constants 
of the large amount of data in this field [8] oer 
2 (cos Gevenbauer polynomial 
indicates a considerable spread, with values of ry : 
D — drag of spher 


Fa) 

treatment. It appears that in range 0-05 to constant functions 

0-40. the value of K may not be uniquely deter- ' j 


mined by ¢ alone, but that particle arrangement ondies 


K up to 100 per cent lower than predicted by our 


modified Bessel function, first kind, zero 


may also be an important variable. modified Bessel function, first kind, first 
The present result Is considere d to be of inter« st order 


since it involves an exact solution of the boundary correction factors to drag as determined 
value probl m proposed by various authors as a by Stokes’ Law (equation 22 and equation 


model of a sedimenting assemblage. It appears to 
odities bess ctio second k 
he in good agreement with analogous Bessel function md kind 
order 
treatments and throws light on the extent of 


disturbance of fluid motion near particles in 


sedimenting systems. It is demonstrated that im . inteser 


cell models for predicting sedimentation dynamics P. (eos 0) Legendre polynomial 


the shape of the outer fluid boundary assumed spherical co-ordinates 

is not important up to substantial solids radius of sphere 

concentration. ae integrals defined by equations (13), (14) 
uniform velocity 
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Letters to the Editors 


Note in connexion with the paper *‘ Kinetics of watergas conversion reaction ’"* 


In his discussion on the possibility of reduction of the 
catalyst (p. 137) the author calculated the oxidizing 
potential for the reaction 5FeO } O, = FeO, at 500 ¢ 
From the thermodynamic point of view, however, FeO 
is not stable below 565 [1]. The reduction of Fe,0, should 
therefore proceed to Fe. If the author has extrapolated 
the iron-oxygen data to below the quadruple point of 
565 C, the oxygen pressure for the FeO — Fe equilibrium 
should have been higher than that for the Fe,0, FeO 
equilibrium. llowever in Fig. 1 the two equilibrium 
pressures are in the natural order as would be the case 
only above the quadruple point. It is concluded therefore 
that the values for the two oxygen pressures as indicated 
in Fig. 1 cannot be both correct. Since the author's 
references to the sources of his thermodynamic data are 
incomplete, it is impossible to verify his calculations 


If a mixture of Fe nap and Fe is contacted with CO and 


H,O under conditions where equilibrium is obtained, at 


* P. Borrouint Chem. Engng. Soc, 1958 9 135. 


S0OC no conversion of Fe,O, in Fe or vice versa should 
occur for a H,O CO ratio of 0-614, using data by 
EmMMeETY and S« Ht The corresponding saturator 
temperature is 75°C, 

The author's experimental value for the ratio at which 
ho reaction occurs ts 0-76, corresponding to 78 C saturator 
temperature Although the author interpreted this 
experiment as a determination of the Fe 30, — FeO equili- 


brium, it must have been the FeO, — Fe equilibrium. The 


discrepancy in the H,O CO m4, might be caused by 
experimental errors (adjustment of saturator temperature, 
measurement of catalyst temperature) or by a diffusional 
resistance. Hlowever if these arguments do not offer a 
satisfactory explanation, the hypothesis that equilibrium 
is established, on which the author's further considerations 
are based, might be at fault. 

G. VAN DEN Bere 
N. V. Mekog, Ijmuiden, Netherlands 
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Tue scope of our work was that of examining if th 
mechanism determining the kinetics of water-gas conver- 
sion was that of the diffusion of the reagents and products 
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numerical calculation is right and, if the results differ from 
those given by Emmet and Scuuvrz, the fault is due to 
the inaccuracy of the thermodynamic data used. The 
casual coincidence of the conditions, thermodyvnemically 
foreseen by us, for the transformation Fe,O, FeO. with 
the conditions experimentally detected for what Dr. van 
DEN Bere thinks rightly to be the Fe transforma- 
tion, has led us astray. We mention once more that 
we are only interested to settle the limit below which the 
value of the ratio H,O CO should not go. What happens 
below that limit is of no interest to us at all and thus we 
have not bothered to check the constitution of the catalyst 
under conditions different from the experimental ones. 
We wish, however, to kindly thank Dr. vAN DEN Bere for 
having clarified our misunderstanding 

P. Borrouint 
Via Tolentino 19, Milan. Italy. 
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